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Averaged controllability for Random 
Evolution Partial Differential Equations 

Qi Lii* and Enrique Zuazua^ 


Abstract 

We analyze the averaged controllability properties of random evolution Partial Dif¬ 
ferential Equations. 

We mainly consider heat and Schrodinger equations with random parameters, al¬ 
though the problem is also formulated in an abstract frame. 

We show that the averages of parabolic equations lead to parabolic-like dynam¬ 
ics that enjoy the null-controllability properties of solutions of heat equations in an 
arbitrarily short time and from arbitrary measurable sets of positive measure. 

In the case of Schrodinger equations we show that, depending on the probability 
density governing the random parameter, the average may behave either as a conser¬ 
vative or a parabolic-like evolution, leading to controllability properties, in average, of 
very different kind. 


Key Words, random evolution Partial Differential Equations, averaged controllability, 
averaged observability, Schrodinger equation, heat equation. 


1 Introduction 

We analyze the problem of controlling systems with randomly depending coefficients in 
the context of evolution Partial Differential Equations (PDEs). More precisely, we consider 
the problem of averaged controllability which consists, roughly, of controlling the averaged 
dynamics, with respect to the random parameters. This problem was introduced and solved 
in [571 in the context of finite dimensional systems, where the same issue was also formulated 
for PDE. 

When the dynamics of the state is governed by a pair of random operators (determining 
the free dynamics and the control operators, respectively), generally speaking, controlling 
the system would require to know the actual value of the random parameters. But this is 
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unfeasible. To avoid this paradox, controls should be independent of the unknown param¬ 
eters. But, of course, this restricts our ability to deal with the randomness of the system. 
Accordingly the control requirement needs to be relaxed. The most natural relaxation is 
to require the control to perform optimally in an averaged sense. This amounts to control¬ 
ling the mathematical expectation of the solutions, making a robust compromise of all the 
possible realizations of the system for the various possible values of the random parameters. 

In this work we consider mainly the heat and Schrodinger equations, with diffusivity and 
dispersivity operators depending on a random variable in a multiplicative manner. We show 
that, while the average of heat equations leads also to a heat-like dynamics, the behavior 
of the averages for the Schrodinger equations depends in a very sensitive manner on the 
density of probability of the random variable so that, in some cases, by averaging, this leads 
to a dynamics of conservative nature, similar to the original Schrodinger equation under 
consideration and, in others, to a parabolic-like behavior. 

Our method of proof combines using Fourier decomposition methods to identify the 
averaged dynamics to, later, utilizing the existing tools developed for the controllability of 
parabolic and conservative systems, to deduce the averaged controllability results. When the 
resulting averaged dynamics is of parabolic nature, averaged null controllability is proved 
for arbitrarily short time intervals and from measurable sets of positive measure. On the 
contrary, when the averaged dynamics is of conservative type, averaged exact controllability 
is proved under suitable geometric conditions on the support of the controls that are by now 
well known in the context of wave-like and Schrodinger-like equations. 

In order to illustrate the effect of averaging and how it may change the dynamics of the 
original system, let us consider the simplest transport equation 

y t + a ■ X 7 y = 0 in R d x [0, oo), 

2/(0) = 2 /o in R d . 

Here 2/0 £ L 2 (R d ) and a(-) : —>■ R d is a d-dimensional standard normally distributed 

random variable, with the probability density 

1 M 2 

p(a) = - je 2“ for a G R d . 

( 27 t) 2 

The solution to ( 1 . 1 ) reads 

2/(x, t, oj ; 2/0) = Vo{x — ta ) for (x, t) G R d x [0, 00). 


Then, the mathematical expectation or averaged state 


y(x,t) = [ y(x,t,u;yo)dF(uj) — ———j f 
Jn (27r) 2 jRd 

(*-£ 




2/o(x — at)e 2 da 


( 27 t) 2 t d Jm d 

solves the following heat equation 


2 /o (z)e 21 2 dz 


Vt ~ jA y = 0 in R d x [ 0 , 00), 

2/(0) = 2/0 hr M d , 


( 1 . 2 ) 
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namely, (p(x,t) = y(x, y / 2 1) solves 


j (p t — Aip = 0 in x [0, oo), 

\ 99 ( 0 ) = y 0 in R d . 

The fact that averages of transport-like equations may enjoy enhanced regularity prop¬ 
erties, first discovered in [1], is well known in different contexts. There are some different 
presentations of this smoothing properties, referred to as “averaging lemmas” in the context 
of kinetic equations (see [9, 21, 22] for example). These smoothing properties are useful 
when proving existence of solutions of linear and nonlinear kinetic equations [6, 12, 13, 39]. 
The proofs of these “averaging lemma” usually employ the Fourier transform in all space- 
velocity-time variables. In [9], the authors gave a proof that only uses the Fourier transform 
with respect to x and v, thus leading to a Fourier representation of the average, evolving in 
time, which is very similar in spirit to the method we employ in this paper. Note however 
that, for control purposes, we need to identify the nature of the evolution associated with 
the averaged variable, and not only its smoothing properties. 

Remark 1.1 The above computation shows that one can get diffusion processes by averaging 
a simple random convection process with respect to its velocity. This is also well known from 
a different perspective, in the context of chaotic and stiff oscillatory systems, that can be 
regarded as the characteristic systems of transport equations (see [11, 19, 45]). 

This example shows that, by averaging, the solution of a random transport equation may 
lead to a solution of a heat-like equation and, consequently, that time-reversible systems may 
become strongly irreversible through averaging. Furthermore, this occurs with the normally 
distributed random variable which is ubiquitous in nature due to the central limit theorem, 
which states that the mean of many independent random variables drawn from the same 
distribution is distributed approximately normally, irrespective of the form of the original 
distribution. Accordingly, in the real world, physical quantities that are expected to be the 
sum of many independent variables, such as measurement errors, often have a distribution 
very close to the normal distribution (see [10]). 

The “smoothing by averaging” effect mentioned above has important consequences from 
a control theoretical point of view as well. Indeed, while the transport equation (1.1), 
for a given value of the random variable a (which determines the velocity of propagation 
of waves), enjoys the property of exact controllability in finite time, proportional to the 
travel time of characteristics to get to the control set (the boundary or an open subset of 
the domain), the averaged heat dynamics is controllable to zero (or any other sufficiently 
smooth target) in an arbitrarily short time and from any subset of the domain where the 
dynamics evolves, without any geometric condition on the support of the controls, involving 
the propagation of characteristics. Accordingly, through averaging, we encounter on a single 
model, with randomly depending coefficients, the classical dichotomy arising in the context 
of controllability of hyperbolic versus parabolic systems (see [56]). 

This paper is devoted to systematically addressing these questions in the context of heat 
and Schrodinger equations. Our aim here is not, by any means, to systematically address 
all the possible scenarios but simply to highlight some of the most fundamental phenomena 
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illustrating how, the existing tools for the analysis of the controllability of PDE, can be 
employed in this averaged context too. It is important to highlight, however, that the 
averaged states do not obey a PDE, not even a semigroup. The dynamics can however 
be represented in such a way that its main controllability properties can be identified, by 
analogy, with some of the main well-known models, and analyzed by similar techniques. 

In particular, we shall show that the averages of heat equations lead to heat-like dynamics 
that enjoy the null-controllability properties of solutions of heat equations in an arbitrar¬ 
ily short time and from arbitrary measurable sets of positive measure. In order to prove 
these results we employ classical techniques based on Carleman inequalities and the Fourier 
expansion of solutions on the basis of the eigenfunctions of the Laplacian generating the 
dynamics. 

In the case of Schrodinger equations we show that, depending on the probability density, 
the average may behave either as a conservative or a heat-like evolution, leading to con¬ 
trollability properties, in average, of very different kind. When the obtained average is of 
parabolic nature the techniques above, employed to treat the control properties of parabolic 
averages, can be applied. However, when the average behaves rather in a conservative way 
we employ specific techniques for the control of wave-like equations. 

The paper is organized as follows. In section 2 we present all these problems in an 
abstract setting in which different relevant PDE models enter naturally. In Section 3, we 
study the null and approximate averaged controllability problems for a class of random heat 
equations. In Section 4, we study the null and exact averaged controllability problem for 
a class of random Schrodinger equations. In Section 5 we give some further comments and 
open problems. 

2 An abstract setting 

Let T > 0 and P C [0,T] be a Lebesgue measurable set with positive Lebesgue measure. 
Let H and U be two Hilbert spaces. Let V C H be a Hilbert space which is dense in H. 
Denote by V' the dual space of V with respect to the pivot space H. Let (Q, T ', P) be a 
probability space. Let {^(a;)}^^ be a family of linear operators satisfying the following 
conditions: 

1. 4(-)ei 2 (fi;W),E)); 

2. A(u>) : D(A) —> H generates a Co-semigroup {S{t,u>)} t >o on both H and V for all 
uj E 

3. S(t, ■)y E L 1 (J2; V) for all y 6 V and t E [0, T\. 

Let B(-) E L 2 (yt-C{U,V)). 

Consider the following linear control system 

f Vt(t) = A(u)y(t) + XE(t)B(u)u(t) in (0, T], 

1 1/(0) = Vo, 

where y 0 E V and u(-) E L 2 (E ; U) is the control. 
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( 2 . 1 ) 


In what follows, we denote by yf,u>-,yo) the solution to (2.1), which is the state of the 
system. Although the initial datum y 0 E V and the control u(-) are independent of the 
sample point u, the state y(t,uj;yo) of the system depends on uj nonlincarly. 

According to the setting above, for a.e. uj E hi, there is a solution y(-,u] yo) E C([0, T]\ V) 
and the expectation or averaged state f Q y(-, uj; y 0 )d¥(uj) E C([0,T]; V). 

We introduce the following notions of averaged controllability for the system (2.1): 

Definition 2.1 System (2.1) is said to fulfill the property of exact averaged controllability 
or to be exactly controllable in average in E with control cost C > 0 if given any yo,yi G V, 
there exists a control u(-) E L 2 (E ; U) such that 

\ u \l 2 {e-,u) < C{\yo\v + \yi\v) ( 2 . 2 ) 

and the average of solutions to (2.1) satisfies 

[ y(T,u>] y 0 )dP(ui) = y\. (2.3) 

Jn 

Remark 2.1 The notion of exact averaged controllability was first introduced in [57]. A 
full characterization was also given in the finite-dimensional setting. In [40], this issue was 
discussed for systems involving finitely-many linear parametric wave equations. 


Definition 2.2 System (2.1) fulfills the property of null averaged controllability or is null 
controllable in average in E with control cost C if given any initial datum y 0 E V, there 
exists a control u E L 2 (E;U) such that 

\u\l 2 (e-,u) < C\y 0 \ v (2.4) 

and the average of the solutions to (2.1) satisfies 

[ y(T,u}-,y 0 )d¥(uj) = 0. (2.5) 

Jn 

Definition 2.3 System (2.1) fulfills the property of approximate averaged controllability or 
is approximately controllable in average in E if given any yo,yi € V and e > 0, there exists 
a control u £ E L 2 (E ; U) such that the average of solutions to (2.1) satisfies 


y{T,iw,y 0 )dF(uj) - y x 


< £. 
V 


Remark 2.2 As in the finite dimensional context ([57]), we can also consider the averaged 
control problem with random initial data, i.e., yo E L 2 (Q]V). Nevertheless, according to 
Remark A.l, this does not lead to any essential new difficulty. Thus, for the sake of simplicity 
of the presentation, we only deal with the case where yo is independent of uj. 


These notions are motivated by the problem of controlling a dynamics governed by a pair 
of random operators (A(u;), B(oj)), where the effective value of the parameter uj is unknown. 
Then, one aims at choosing a control, independent of the unknown uj, to act optimally in 
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an averaged sense, making a robust compromise of all the possible realizations of the system 
for the various possible values of the sample point c o. Similar problems can be considered in 
the case where the initial datum to be controlled depends on c o too. 

We have introduced the notions of exact/null/approximate averaged controllability in the 


framework of random evolution equations but similar concepts make sense for parametrized 
evolution equations (see [57] for example). In that context it is sufficient to replace the 


expectation by a weighted average of the parameter-depending controlled states. 

Remark 2.3 In the present context of randomly depending operators, the classical subordi¬ 


nation properties of some control properties with respect to the others, that are classical for 


a given system, have to be addressed more carefully. Of course, averaged exact controllabil¬ 
ity implies the averaged null and approximate controllability properties as well. But, when 
A and B are independent of co and A generates a Co-group, exact controllability is also a 
consequence of null controllability. However, the later may fail when considering averaged 
controllability properties, as shown in the example in Remark f.6 below. 

Remark 2.4 For parametric control systems one can also consider the problems of simulta¬ 
neous controllability and ensemble controllability, which concern the possibility of controlling 
all states with respect to different parameters simultaneously by one single control. We refer 
the readers to [37, 38] and [8, 35] for an introduction to these notions, respectively. Of course 
the property of averaged controllability we consider here is weaker than these other ones since 
we only deal with the average of the states with respect to the parameters. But, as we shall 
see, averaged controllability properties may be achieved in situations where simultaneous and 
ensemble controllability are impossible. For instance, let us consider the following control 
system: 



( 2 . 6 ) 


where 



By Theorem 1 in [57], we know that the system (2.6) is null averaged controllable. 


But system (2.6) is not simultaneous null controllable. Otherwise there would exist a 
u G L 2 (0,T) such that 







which implies that yo — 0. 

Remark 2.5 The connection between averaged and simultaneous controllability was analyzed 
in [fO] through the following optimal control problem: 


Minimize 
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for all 


u E lu E L 2 (0,T;U) : The corresponding solution y satisfies that / y(T, co)dF(co) = yi 


In that paper it is proved that for every k, J k (-) has a unique minimizer u k (-) and that, 
if H and U are finite dimensional and the system (2.1) is simultaneously controllable, then 
{u k }^ =1 weakly converges to a simultaneous control. 

Remark 2.6 In the particular case that A is independent of to, the averaged controllability 
problems can be reduced to the classical controllability ones by setting 

B = [ B(u)dF(u), y(t )= / y(t,u;y 0 )dP(u). 

Jn Jn 

Then we have that 

f Vt(t) = Ay(t) + XE(t)Bu(t ) in ( 0 , T ], . , 

\ 27(0) = y 0 . 

The exact (resp. null, approximate) averaged controllability problems of (2.1) are equivalent 
to the exact (resp. null, approximate) controllability problem of (2.7). 

Remark 2.7 Random evolution equations can be used to model lots of uncertain physical 
processes (see [49, 50, 51] for example). Several notions of controllability have been introduced 
but, as far as we know, all of them concern driving the state to a given destination by 
a control depending on co (see [27, 43, 44] an d the references therein). The property of 
averaged controllability is, however, independent of the specific realization of tv. 


Following the classical approach to deal with controllability problems, we introduce the 
adjoint system, which also depends on the parameter co: 

I ~zt(t) = A*(u)z(t) in [0, T), (2 

1 z(T) =z 0 , [ ' J 


where z 0 E V'. 

Note that, in this adjoint system the initial value (at time t = T) is taken to be inde¬ 
pendent of to. This is due to the fact that, although y(T,uo\yfi) depends on co, its average, 
which belongs to V, is of course independent of co. Then, to deal with the averaged state 
it is sufficient to use as test functions, adjoint states departing from configurations that are 
independent of co. This is the reason we choose the final datum of (2.8) to be independent 
of co. 

As dual notions of the properties of averaged controllability above we introduce the 
following three concepts of averaged observability. 


Definition 2.4 System (2.8) is exactly averaged observable or exactly observable in average 
in E if there is a constant C > 0 such that for any z 0 E V', 


|-o| 2 v<C XE(t) / B(uj)*z(t,uj]Zo)dF(uj) 


dt. 


u 


(2.9) 
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Definition 2.5 System (2.8) is null averaged observable or null observable in average in E 
if there is a constant C > 0 such that for any z 0 £ H, 


z(0, to; £o)dP(a;) <C / Xe{1) / B(uj)*z(t, to; z 0 )dP(uj) 

• V' In .In 


dt. 


u 


( 2 . 10 ) 


Definition 2.6 System (2.8) is said to satisfy the averaged unique continuation property in 
E if the fact that xe f n B(oj)*z(-,uj] zo)dP(u>) = 0 in L 2 (0,T;t/) implies thatzo = 0. 


The adjoint system (2.8) being exactly ( resp. null) averaged observable means that one 
can estimate the norm of the average of the final (resp. initial) data of the adjoint system, 
out of partial measurements done on the averages of the adjoint states, with respect to oj. 
These concepts have their own interest when dealing with the observation of random systems. 
The actual realization of the system depending on uj being unknown, it is natural to address 
the problem based on the measurements done on averages. 

The weakest notion of averaged observability under consideration is averaged unique 
continuation. System (2.8) satisfies the averaged unique continuation property when its 
state can be uniquely determined by the partial measurements done on the mathematical 
expectation. It is a natural generalization of the unique continuation property of evolution 
equations. 

The average of the adjoint state, being represented by f Q z(t, u; zo)dP(u), does not satisfy 
the semigroup property and it is not a solution to an evolution equation. Thus, one can not 
directly employ the existing results on the observability of PDEs to establish the averaged 
observability of the adjoint system (2.9). However, as we shall see, by carefully analyzing 
and identifying the dynamics generated by the averages of the adjoint states, we shall be 
able to apply the existing PDE techniques. 

In this paper, we mainly consider the case that A(u>) = a(uj)A and B(u>) = B, where A 
generates a Co-semigroup on H , a(co) is a random variable and B £ C(U, H ). We will show 
that the controllability properties of the system (2.1) depend on the choice of the random 
variable. We only consider the following commonly used ones: 

1. Uniformly distributed random variable , with probability density function p(-) on [a, b], 
where 0 < a < b, a(-), given by 


p(a) 


1 

b — a’ 

0, 


if a £ [a, b], 

if a £ (—oo, a) U ( b , oo). 


2 . 


by 


Exponentially distributed random variable , with probability density function p(-) given 


p(a) 


e c \ if a > c, 
0, if a < c, 


( 2 . 11 ) 


for a given positive number c. In what follows, for simplicity, we choose c = 1 in (2.11). 

3. Standard normally distributed random variable , with probability density function of 
p(-) reads 








4. A random variable with standard Laplace distribution , with density function p(-) given 
by 

p(a ) = for a G M. 

5. A random variable with standard Chi-squared distribution , with density function p(-) 
given by 


p(or,k) = < 

a 2 1 e 2 

O 

A 

3 

2 ir(f) ’ 


\ 

l o, 

if a < 0, 

5 — — 1 — Ol d 

az e aa. 



6. A random variable with standard Cauchy distribution , with density function /?(•) given 


by 


p(a) = 


1 


7r(l + a 2 ) 


for a 6 


3 Averaged controllability for random heat equations 

In this section, we study the null averaged controllability problem for random heat equa¬ 
tions. Of course, the property of averaged exact controllability has to be excluded because 
of the very strong regularizing effect of heat equations, that is preserved by averaging. Thus, 
we focus on the properties of the null and approximate averaged controllability. 

Let T > 0 and G C (d E N) be a bounded domain with C 2 boundary dG and consider 
the following controlled random heat equation: 


{ y t - aAy = Xg 0 xeu in G x (0, T), 

y = 0 on 9G x (0,T), (3.1) 

?/(0) = y 0 in G. 

Here yo E L 2 (G), and Go and E are subsets of G and [0, T] respectively, where the controls 
are being applied. The constant diffusivity a : H —» M + is assumed to be a random variable. 
In this section, we make the following assumptions on Go and E: 

(Al) G 0 C G is a nonempty open subset. 

(A2) E C [0, T] is a Lebesgue measurable set with positive measure. 

We first analyze the dynamics of the mathematical expectation (average) of the solutions 
of the random heat equations under consideration. This will be done using Fourier series 
expansions and will give us an intuition for the averaged observability and controllability 
results to be expected and that will be proved below. 

Consider first the following uncontrolled random heat equation: 

[ Vt~ aAy = 0 in G x (0, T), 

\ y = 0 on dG x (0, T), (3.2) 

I 2/(0) = llo in G. 
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We decompose solutions in Fourier series on the basis of the eigenfunctions of the Dirichlet 
laplacian. To be more precise, consider the unbounded linear operator A a on L 2 {G) given 
as 

f D(A A ) = H 2 (G)nH 1 0 (G ), 

I A a u = —Ait, for any u £ D(A a ). 

Let us denote by {A J }°L 1 (with 0 < Ai < A 2 < • • •) the eigenvalues of A a and let be 

the corresponding eigenfunctions such that \ej\ L 2 ^ = 1 for j £ N. 

We assume that the initial datum takes the form y 0 = Vo,j e j *= L 2 (G). The averaged 
state can also be described in Fourier series as follows, distinguishing the probability densities 
above. 

Case 1. If «(•) is a uniformly distributed random variable on [a, b], where 0 < a < b, 
then, 



Remark 3.1 The values of a uniformly distributd random variable, which is a relevant one 
in 'practice, are uniformly distributed over an interval, i.e., all points in the interval are 
equally likely. It models the random phenomenon with “equally possible outcomes”. When 
there is no any a priori knowledge for a other than a(u) £ [a, b] for all uj £ 12, this is the 
best possible choice. More details can be found in [26], 


Case 2. When a(-) is the exponentially distributed random variable, 


y(x,t,uj-,y 0 )d¥(uj ) 



y^/!Jo,je Xj0lt ejda 
3 =1 


oo 


= £ 


Xjt 


+ 1 


Vo,j e Xjt ej. 


(3.4) 


Remark 3.2 The exponentially distributed random variable is one of the most important 
random variables that can be used to describe the time between events in a process where 
they occur continuously and independently at a constant average rate. 


In both cases, the mathematical expectation of the solution of the parameter-depending 
heat equation evolves according to a heat-like dynamics. The representation of the averages 
on the basis of eigenfunctions exhibits the exponential decay and smoothing effects that are 
prototypical of heat-like problems. 

Accordingly, the following null averaged controllability result holds. 

Theorem 3.1 Let (Al) and (A2) hold. Assume that, either «(•) is a uniformly or expo¬ 
nentially distributed random variable. Then, the system (3.1) is null controllable in average 
with control u(-) £ L 2 (0,T; L 2 (G 0 )). Further, there is a constant C > 0 such that 

\u\l 2 {o,t-,l 2 (g 0 )) A C||/o|l2(g)- (3-5) 
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Remark 3.3 In the context of heat equations, the null controllability result with controls 
supported in measurable sets is related to the bang-bang property of the time optimal control 
problems with constrained controls (see [48, 52] for example). 

More precisely, let us consider the following heat equation: 

{ Vt- Ay = xg 0 u in G x [0, +oo), 

y — 0 on dG x [0, +oo), (3.6) 

2/(0) = Vo in G. 

Here the initial state f/ 0 G L 2 (G) and the control is assumed to belong to the constrained set 
of admissible controls 

u G Um = {«G L°°(0, +oo; L 2 (G )) : |m|l 2 (G) < M a.e. t G [0, +oo)} 

for some M > 0. Let T* be the min {t : y(t;u,yo) = 0, u G Um}- A control u* G Um such 
that y(T*-,u,y 0 ) = 0 is called a time optimal control and satisfies the bang-bang property if 
|^(^)|l 2 (G) = M for a.e. t G [0,T*]. The proof of this fact requires the null controllability 
with controls supported in measurable sets (see [52]). 

The same bang-bang problem can be formulated in the context of averaged null controlla¬ 
bility we are considering here. However, the techniques in [52] do not seem to apply because 
the averages of the heat processes under consideration do not satisfy the semigroup property. 


To prove Theorem 3.1, as in the deterministic frame, we introduce the following adjoint 
system: 

{ z t + a A z = 0 in G x (0, T), 

£ = 0 on <9G x (0, T), (3.7) 

z(T) = z 0 in G, 

where Zo G L 2 (G). As mentioned above in the abstract setting, the initial data (at time 
t = T) of the adjoint system are assumed to be independent of the random parameter. 
According to Theorem A.2, we only need to prove that (3.7) is null observable in average, 
which is a corollary of the following result. 

Theorem 3.2 Let (A1 ) and (A2) hold. Assume that «(•) is either an uniformly distributed 
or an exponentially distributed random variable. Then, there exists a constant C > 0 such 
that for any yo G L 2 (G ) it holds that 


z(0,uj-,y o )dF(uj) 


in 


L 2 (G) 


< c 


z(t, yo)dF(uj) 


L 2 (Go) 


dt. 


(3.8) 


An immediate corollary of Theorem 3.2 is as follows: 

Corollary 3.1 Let (Al ) and (A2 ) hold. Assume that a(-) is either an uniformly distributed 
or an exponentially distributed random variable. Then the system (3.7) is null observable in 
average. 
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Proof : From Holder’s inequality, we have that 


z(t, uj; z 0 )dP(uj) 


L 2 (G 0 ) 


dt < 


dt 


< \fm(E ) 


/ z(t,ur, Zo)dF(u>) 

'n 

z(t, oj; zo)dP(w) 


LHG 0 ) 

2 

dt 

L 2 (Go) 


dt 


This, together with the inequality (3.8), implies that the system (3.7) is null observable in 
average. □ 

To prove Theorem 3.2, we adopt the strategy developed in [46] in the context of the heat 
equation, using spectral decompositions. The null averaged observability inequality is built 
in an iterative manner. More precisely, we decompose the set E into an infinite sequence 
of connective (in time) subsets in which an increasing number of Fourier components of the 
average of the solution is observed with uniform observability constants. By iteration, the 
final datum is observed. To apply this strategy, we need to use classical results on how to 
divide E into an infinite sequence of subsets of positive Lebesgue measure. We also need 
to know how to observe a finite number of the Fourier components of the average of the 
solution. These ingredients are given in the following two lemmas. 


Lemma 3.1 [48, Proposition 2.1] Let E C [0, T] be a measurable set of positive Lebesgue 
measure m(E). Let £ be a density point of E. Then for each a > 1, there exists an £ i E (£,T) 
such that the sequence {4}*4 1; given by 


~k +1 — + 


h-£ 


(3.9) 


satisfies 


m(E fl (4 + i,4)) > 


-k — <T+1 


(3.10) 


Lemma 3.2 [fl, Theorem 1.2] There is a constant C\ > 0 such that for any r > 0 and 

{dj}\j<r C 

2 


]T \ aj \ 2 j 2 < c ie Ci ^( / ( 1 eA x ) 


Go 


dx 


(3,11) 


We are now in conditions to proceed to the proof of Theorem 3.2. 

Proof of Theorem 3.2: We only give a proof for the case where a(-) is an exponentially 
distributed random variable. The proof for that «(•) is a uniformly distributed random 
variable is very similar. 

Put z(x,t ) = f n z(x,T — t,u] zo)dF(cu). Let zoj = {zq, ej)ip-(G)- Then, similar to the 
computation of (3.4), we have that 

OO ^ 

z(x,t) = ^ X t + 1 Zo h e ~ Aj V 
j =i J 

We only need to prove that 


OO 

Y ' V7f ./ 

3 = 1 

< c 

L 2 (G) 


12 


z(x,t,ui)\v,( 0a )dt. 


(3.12) 














Note that, in the right hand side of this inequality, the observation is done in the L 1 ( E ; L 2 (Gq))- 
norrn. Thus, the result is even stronger than the one we actually need, in which the obser¬ 
vation is done in L 2 (E; L 2 (G 0 )). As a consequence of the inequality above we shall prove, 
the controls we shall build will belong to L°°(E ; L 2 (Gq)). 

Let X r = span for each r > 0. For any £ G L 2 (G), we put 

OO 

s ( i ’fl = Ev7TT^ A '‘ ei ’ (3 ' 13) 

j =1 2 

where £j = (£,&,■) l 2 (G)- Then we see that 

|^(^0|l 2 (g) < |^(s,0|l 2 (g) , for 0 < s < f < T, (3.14) 


and 


S((,0|i S(O ) < e r< * S ^|5(s, 0 Il 2 (G). for all Z £ V. and 0 <s<t<T. 


(3.15) 


Let £ be a density point for E. By Lemma 3.1, for a given a > 1, there exists a sequence 
{£k}kLi satisfying (3.9) and (3.10). 

We now define a sequence of subsets {E k }™ =1 of (0, T) in the following way: 


Eu = It — 


ofc + l 


: t G E fl ( l k - i-i + 


-k +1 


, 4 ) f, for k G N. 


Clearly, E k C (4+i, 4+i + | (4 — 4+i))- From (3.10), we have that 


m(-Efc) = m[E fl (4+i + 


-k +1 


6 


= m Ffl 


(4+i) 4) \ K fc+i) 4+i + 


4 4+i \ 

— 6 y 


> m(A n (4+i,4)) 
4 — 4+i 


tfc+i 


> 


6 


Let b > a be a positive number such that 


b C\ + 61n(12Cia) 

a (a — 1)(4 — 4 


Set r k = & 2fc . From (3.14), we have that for any £ G X rk , 

/■^*;+i+|Ffc—4+i) 

J £k +1 


Xfi fc (4 5^4+i + -(4 - 4+i),£) 


L2(G) 


dtt 


f-^fc+i+f ^fc+i) 


< 


'h+i 




(3.16) 


(3.17) 


(3.18) 


13 














From (3.11) and (3.17), we find that 


4 4+i 


S (4+1 + g(4 — 4+04) 
< m(E k ) s(i k +i + jj(4 - 4+i)4) 


L2(G) 
L2(G) 


< 


f4+l+§(hfe — 4+l) 


' 4+1 


XE k (.t)\S(t,£)\ L2{G) dt 


(3.19) 


< C x e Cl ^ 


f4+l + f(4 — 4+l) 


' 4+1 


XB fc W|5'(t4)| i 2 (Go) ^- 


Let z 0 = Zq + Zq, where z\ £ X rk and z% £ X^. Taking £, = S k g fc+1 , ^) hr (3.19), we get 
that 


— <++1 


6 

r4+i+§(4—4+i) 


|S(4,^o)U 2 (g) 


< 


* tk+i 


XE k {t ) 


5(t+^—^±1,4 




L 2 (G) 


dt 


< C x e Cl ^ 

< C x e Cl ^ 


r4+i+§ (^fc—■4+i) 


' 4+1 

r^k 

'4+1“ 


X^(0 


St + 


4; 4+1 1 


6 


> 4 


L 2 {Go) 


dt 


4 "4+1 


U - 


<-fc — <-fc+i 


4(L-o)|l2(Go)^- 


By (3.16), we have that 

4 4+i 

6 


XEfc t ~ 


= XE(t ), for any t £ ( 4+1 + 


4 4+i 

~6 


Combining (3.14), (3.15), (3.20) and (3.21), we find that 


4 4+i 

6 

< C x e Cl ^ 

< Cie Cl ^ 

< C x e Cl ^ 


|*S'(4,^o)U 2 (g) 


r4 


L , 4“4 +1 
'4+lH-g- 


/■4 

'4+i 

r4 


4~4+i 


h k+1+ l -^±L 


XE(t)\S(t,zl)\ L 2 {Go) dt 

XE(t)(\S(t,Z 0 )\ L 2( Go) + \S(t,Z%)\ L 2 { G))dt 
XE(t)\S(t, Z 0 )\ L 2( Go )dt 


+C ie Cl ^( 4 - 4 + 1 ) 


o l n ,4 4+1 2 

*j ( 4+i H---, ^0 


L2(G) 


r4 


<C ie Gl ^ XE(t)\S(t,z 0 )\ L 2 {Go )dt 
J 4 + 1 

+C x e c ^{ 4 - 4+i)e^^ ±irfc |^(4+i^ 0 2 )| L 2 (G) - 


(3.20) 


(3.21) 


(3.22) 
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Therefore, we obtain that 


4 6 <t+1 |S(4,3))k»(G) 

<Cye c [\ E (t)\S(t,z 0 )\ LHCo) dt+C 1 e c ^(e i ~£ M )e-' J ^^\S(e M ,zi)\ L , (aj 

J £k+l 

+ 4 ~ 4+1 |S(4,^)U» (C ) 

f XE(t)\S(t,z 0 )\ LHCo) dt+C\e c ^(e,-t M )e-‘-^'‘\S({ M ,4)\ LHG) 

■'4+1 

+ 4 ~4h- 1 e - rt( 4-4«) |g ( 4+1 ’ z g)| L , (c) . 

(3.23) 


Thus, it holds that 


4 — -tfc+i 
6 


|S(4, ^o)|l 2 (g) 


r4 


<C x e Cl ^ XE(t)\S(t,z 0 )\ L 2 { G o) dt 

d (-k+i 

+ (4 - 4+i)e- f ^ ±irfc (4 1 e Cl ^ + 1) |S'(4+i, *o)| 


L 2 {G)' 


This concludes that 

4 — 4+i | Q ,„ \ i Cie 01 ^ -f i _ 4-4+i 

^=|S(4,-o)|l 2 (g) - Cie G lv /rT (4-4+i)e 6 |S(4+1, ^o)U 2 (G) 

/■4 

< / XE(t) |S(t, £ 0 ) |l 2 (G 0 )^- 

J 4+1 

By summing the inequality (3.25) from k — 1 to k — oo, we obtain that 


4-4 


|S(4, ^o)|l 2 (G) + ^ 4|S(4+1, ^o)|l 2 (G) < / XE(t)|S(t,2o)| L 2(Go)dt, 

k=1 

fc+i — 4+2 4' 1 e c ’ 1 \ /F fc + 1 


6C'ie c ' 1 v /F >" 
where 

^ k QC\e Cl X rk + 1 

From (3.18) and r*, = b 2k , we have that 


-(4 - 4+i)e _jLZ 6 ttirfc , 4 = 1,2,--- 


(3.24) 


(3.25) 


(3.26) 


fk > 0 for any k = 1,2, -- 
This, together with (3.26), deduces that 

|S(4, 2o)|l 2 (G) < —d—j—J E \S(t, z o)\ L 2^ Go) dt. 

Since 4 < T, we can find a constant 4 > 0 such that 

C 

- -—— > e~ Aj ( T -4 for every j G N. 

1 + Ajti 


(3.27) 


(3.28) 
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From (3.27) and (3.28), we get that 


J2 z °’i e 

3 = 1 


L 2 (G) 


<c[ \S(t,Zo)\ L 2(Go) dt - 

J E 


This completes the proof. 

As an easy corollary of Theorem 3.2, we have the following result. 


(3.29) 

□ 


Theorem 3.3 Let (AlJ and (A2) hold. System (3.1) is approximately controllable in av¬ 
erage, provided that a(-) is a uniformly distributed or an exponentially distributed random 
variable. 


Proof: According to Theorem A.3, we only need to prove that the unique continuation 
property in E is satisfied. 

Assume that z — 0 in Go x E. From Theorem 3.2, we obtain that 


2 ( 0 , 0 ;; 2 0 )dP(o;) 


L 2 (G) 


< C 


z(t , u; 2o)dP(n;) 


L 2 (G 0 ) 


dt. 


Hence, f n z(0, uj; zo)dP(cv) = 0. On the other hand, if zqj = f G zoGjdx, 


n pOO OO OO h 

0 = / 2 ( 0 , uj] z 0 )dF(uj) = / e -0-b ^2 Zoje-^^ejda = ^ — 1 z 0 , j e~ XjT e j . (3.30) 

Jn J 1 j=i j =1 a 3 ' 

Since {ej}ff 1 is an orthonormal basis of L' 2 (G), it follows that x .^ +1 z 0 ,j e ~ XjT = 0 for all 
j G N. Thus, we get that Zqj = 0 for all j G N, which implies that zq = 0. □ 


4 Averaged controllability for random Schrodinger equa¬ 
tions 

4.1 Preliminaries 

In this section, we study the null and exact averaged controllability problems for a class 
of random Schrodinger equations of the form 

{ y t - iaAy = Xg 0 xeu in G x (0, T), 

y = 0 on dG x (0, T), (4.1) 

2/(0) = 2/0 bi G. 

Here a(-) : H —y M is a random variable, the initial datum y 0 belongs to L 2 (G ) and Go is a 
suitable subdomain of G. 

In this time-reversible setting the Schrodinger equation is well-posed whatever the sign 
of a is, contrary to the heat equation. Thus, we have more choices for the random variable 
«(•)• 
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The average of the solutions to random Schrodinger equations may lead to very different 
dynamics, depending on the random variable under consideration. To see this, let us first 
consider the Schrodinger system in the absence of control: 

{ y t — iaAy = 0 in G x (0, T), 

y = 0 on 3G x (0,T), (4.2) 

m = ilo ill G. 

Here y 0 = Vo,j e ] e L 2 (G). 

Case 1. When a(-) is a uniformly distributed random variable on [a, b\, where a, b G M, 
then, 



Case 2. When a(-) is an exponentially distributed random variable, then 


y(x,t,uj-,y 0 )dF(uj) = e 


J 1 

OO 


oo 00 

-(«-!) ^2 yoje-^ejda 
3 = 1 


(4.4) 


Ej 

3 = 1 


■ , — tVo 

i\ i t + l y ’ 3 J 


Case 3. For the normally distributed random variable a(-) we have, 


j y(x,t,u}-,y 0 )d¥(uj) = —= j 


e 2 

°° 3 =1 




(4.5) 


I] yo ; . 

3 = 1 


_1 \2 j2 

W en¬ 


case 4. When a(-) is a random variable with Laplace distribution we have 


y(x,t,uj-,y 0 )dF(uj) 



e~ H yo,je~ iXjat ejda 

3=1 


oo 


= E 


i + 


yo,j e 




3 ‘ 


(4.6) 


Remark 4.1 The Laplace distribution can be thought of as two exponential distributions 
(with an additional location parameter) spliced together back-to-back. It governs the difference 
of two independent identically distributed exponential random variables and can be regarded 
as the generalization of the exponential distribution on the whole real line. More details can 
be found in [29], 
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Case 5. For the Chi-squared distribution it holds: 


y(x,t,uj-,y 0 )dF(uj) = 


/»00 —_ I __ 

f «2 1 e 2 

o 2fr(|) ^ 
l 


E y°’i e lXjate i da 


E- 

(1 + 2i\jt)2 


Vo,je 


—iX-it „ 


(4.7) 


Remark 4.2 A Chi-squared distributed random variable is the sum of the squares of k in¬ 
dependent standard normally distributed random variables. It is one of the most widely used 
probability distributions in inferential statistics. We refer the readers to [25] for more details. 

Case 6 . For the Cauchy distribution we have 


y(x,t,u;y 0 )dP(u) = 


Loo tt(1 + « 2 ) “ 


lXjat Cjda 


= E y°-i e Xjte i- 

j =i 


(4.8) 


Remark 4.3 The Cauchy distribution is associated with many processes, including reso¬ 
nance energy distribution, impact and natural spectral and quadratic stark line broadening. 
It also has important connections with other random variables. For example, when 71 and 
72 are two independent standard normally distributed random variables, then the ratio 71/72 
has the standard Cauchy distribution. More details can be found in [25]. 


Remark 4.4 From (4.8), we know that when a is a random variable with Cauchy distri¬ 
bution, then the average of the solution to the random Schrodinger equation (4.2) becomes 
a solution of the heat equation. This is another example that after averaging, one enjoy 
enhanced regularity properties. 


According to the above results, there are essentially two different dynamics for the aver¬ 
ages, depending on whether the time-exponentials entering in the Fourier expansion are real 
or imaginary. In cases 3 and 6, the average has a heat-like behavior. However, in cases 1, 2, 
4 and 5, the average has a Schrodinger-like behavior. 


4.2 Null averaged controllability 

We first recall the following assumptions on Go and E: 

(Al) Let Go C G be a nonempty open subset. 

(A2) Let E C [0, T] be a Lebesgue measurable set with positive measure. 

Theorem 4.1 Let /Al/ and (A2) hold. If a is a random variable with normal distribution 
or Cauchy distribution, then the system (4.1) is null controllable in average with control 
u(j E L 2 (0, T; L 2 (Go)). Further, there is a constant G > 0 such that 

Ml 2 (o,t ; l 2 (g 0 )) L G|?/o|l 2 (g)- (4.9) 
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Remark 4.5 Note that, in the present case, the random Schrodinger equations is null con¬ 
trollable in average without any assumption on the support G 0 of the control, other than being 
of positive measure. This is in contrast with the well known results on the null controllability 
of Schrodinger equations, where G 0 is assumed, for instance, to fulfill the classical Geometric 
Control Condition (GCC)(see [32, 33] for example) or other geometric restrictions associated 
to multiplier techniques or Carleman inequalities (see [30, 42] for example). In the present 
case, these restrictions on Go are not needed since the averages behave in a parabolic fashion. 

To prove Theorem 4.1, we introduce the adjoint system of (4.1) as follows: 

{ z t + iaAz = 0 in G x (0, T), 

£ = 0 on<9Gx(0,T), (4.10) 

z(T) = z 0 in G. 

By Theorem A.2, we only need to prove the following result. 

Theorem 4.2 Under the assumptions of Theorem 4-1 the system (4.10) is null observable 
in average if a is a random variable with normal or Cauchy distribution. 

Indeed, we have the following stronger observability estimates. 

Proposition 4.1 Let (A.1) and (A2) hold. Assume that «(•) is a random variable with 
normal distribution or Cauchy distribution. Then there exists a constant C > 0 such that 
for any z 0 G L 2 (G), it holds that 


z[f),uj ; 2 0 )dP(a;) 


LUG) 


< C 


z(x, t, u ;; £o)dP(a;) 


LUGo) 


dt. 


(4.11) 


The proof of Proposition 4.1 is very similar to the one for Theorem 3.2. We omit it here. 
Remark 4.6 If a(-) is a random variable with Cauchy distribution, then we have that 

OO 

7^7 z "J e 

(4.12) 


z(x,t,uj] z 0 )dF(uj) = 


00^(1 + a 2 ) “ 




e - A .( T -t) 


e j. 


In this case, Proposition 4-1 is an immediate corollary of the observability estimate for heat 
equations. 

This is an example of a system that is null but not exactly controllable in average. 


4.3 Exact averaged controllability 

In this subsection we consider the cases where the averages behave as Schrodinger-like semi¬ 
groups. We assume that E = [0, T] so that the control is active in any time instant within 
the time interval [0, T]. 
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Let us first consider the following Schrodinger equation 


{ ( p t + KiAtp = 0 in G x (0, T ], 

V = 0 on 9G x (0,T), (4-13) 

¥>(0) = Vo in G, 

where /t£l \ {0} and Vo G L 2 (G). 

We make the following assumption in this subsection on G 0 : 

(A3) Whatever T > 0 and k ^ 0 are, there is a constant C > 0 such that for any 
Vo E L 2 (G), the solution ip(-, •) to (4.13) satisfies 

\Vo\b(G) < C \ I \v\ 2 dxdt. (4.14) 

Jo J Go 

We refer the readers to [3, 24, 32, 47] for the study of (4.14) under different conditions for 
Gq. In those articles one can find various sufficient conditions on the subset Go so that the 
observability inequality above holds, depending of the techniques of proof employed (multi¬ 
pliers, Carleman inequalities, Microlocal analysis). In particular, this observability inequality 
for the Schrodinger equation holds as soon as it is satisfied for the wave equation in some 
time horizon. Thus, in particular, it holds under the classical Geometric Control Condition 
(GCC) guaranteeing, roughly, that all rays of geometric optics enter the observation set Go 
in some uniform time. 

We have the following observability result for the system (4.13). 

Theorem 4.3 The following results hold: 

• Let a(-) be a uniformly distributed random variable on an interval [ a,b\. Then, the 
system (4.1) is exactly controllable in average with V = H = L 2 (G ) and U = H~ 2 (G 0 ). 

• Let a(-) be an exponentially distributed random variable. Then, the system (4.1) is 
exactly controllable in average with H = L 2 (G), V = H 2 (G) fli/g(G) and U = L 2 (G 0 ). 

• Let a(-) be a random variable with Laplace distribution. Then, the system (4.1) is 
exactly controllable in average with H = L 2 (G), V = H 4 (G) D/7g(G) and U = L 2 (G 0 ). 

• Let a(-) be a random variable with standard Chi-squared distribution. Then, the system 
(4.1) is exactly controllable in average with H = L 2 (G), V = H k (G ) ft L7 q(G) and 
U = L 2 (G 0 ). 

Proof of the first conclusion in Theorem 4.3: We divide the proof into two steps. 

Step 1. In this step, we show that the average of solutions can be represented by the 
difference of the solutions of two Schrodinger equations. This was already noticed in [58]. 
We present the argument here for the sake of completeness. 

We let zo G L 2 (G) in (4.10) and z(x,t]Zo) = f Q z(x,T — t,u; zo)dP(u). Assume that 
z 0 = £7= -i z o,j e j- Similar to (4.3), we have that 

1 00 1 

z(x,t;z 0 ) = - - y Zoj ( e -^- e -^ bt ) ej . (4.15) 

b — a ^ iXjt v ' 

3=1 j 
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Let 


E ~o,j 

iX-i 


z o,j -i\. 


at 


3 = 1 


*w>\) = 


E ^o,j 


z 0,j -iXM 


3 = 1 


Clearly, z a (-, •) and z&(-, •) solve the following equations respectively: 

+ aiAz a = 0 in G x [0, T], 
z a = 0 on dG x [0, T], 

z a { 0) = A a 'z q in G, 


Zb,t + biAzb = 0 in G x [0, T], 

Zb = 0 on dG x [0, T], 

z b { 0) = A^z 0 in G. 

tz(-, •) = (z a - z b )(-, •)• 


(4.16) 


(4.17) 


Further, 


Step 2. In this step, we establish the exactly averaged observability estimate. From 
(4.16) and (4.17), we have that 

(id t + aA){id t + bA){z a - z b ) 

= {id t + a A) {id t + bA)z a — (id t + aA){id t + bA)z b 
= {id t + bA){id t + aA)z a = 0. 


Hence, we know that {id t + bA){z a — z b ) solves 

iipt + aA <p = 0 in G x (0, T], 

(/? = 0 on dG x (0, T), 

<^(0) — {b — a)zo in G. 

By assumption (A3), for any zq E L 2 {G), it holds that 


(4.18) 


\ z o\ l 2 (g) A c 
< C 


\<p{x, t)\ 2 dxdt < C 


| {id t + b.A) {z a — Zb) {x, t ) 1 2 dxdt 


'o JG 0 

rT r 


'0 JG 0 

rT 


1 0 •' Go 


\{id t + bA)[tz{x,t] zo)]\ 2 dxdt < C j \tz{x,t; Zo)\ 2 H2 ^ G ^dt (4.19) 


<C \S{x,t] z 0 )\ H 2 {Go) dt. 


The proofs of the second to the fourth conclusion in Theorem 4.3 are very similar. We 
only give that for the second one. 

Proof of the second conclusion in Theorem f.3 : Let z 0 E V' — [H 2 {Q) fliT^O)]' in (4.10) 
and z{x, t; z 0 ) = f Q z{x, T — t,w; z 0 )dF{u>). We only need to prove that 

\ z o\v' <C / \z{x, t) z 0 )\ 2 dxdt. 

JO J Go 
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(4.20) 


We divide the proof into two steps. 

Step 1. In this step, we prove that a “weak” version of the exact averaged observability, 
that is, there is a lower order term in the right hand side of the inequality. 

Assume that z 0 = Yl'jLi Zo,j e j £ V' . Then, 

OO 

z(x,t,cj;z 0 ) = y g 0 je~* aAj ' (r ~* ) e J -. 

5=1 


Similar to (4.4), we have that 


z(x,t;z 0 )= / z(x,T -t,cj;z 0 )dP(uj) = V] ——— -z 0J e lXjt ej. (4.21) 

Jo j=1 lAjt + i 


This implies that for any 5 > 0, 


|z(', S A>)|L 2 (<5,T;L 2 (G)) < C(<J, T)|z 0 |v'. 


(4.22) 


Let 


n(x,t) =^2jy-iZoje * A 'V 


i\d 

5=1 J 


From assumption (A3), for a fixed 5 > 0, we have that 


|2 °L _ I LI ~xd 


5 = 1 


p-^dp . 

J c c 5 


2 <c/ T 


^0,5 —i\-t 


L2 (G) J 5 J Go ' “7 W 


E*0,i 

ir e " v " e 

5=1 


A ~ 


dxdt 


— C I / |ti;(x, t)| 2 da;df. 

J S J Go 


Therefore, 
ko| v 7 


(4.23) 


< C / / |tw(o;, t)| 2 d:r<it 

J<5 ./g 0 

<C / |tA(x, i; + / / |tf(x, f) — tz(x, f; £o)| 2 <icJt 

L Jtf Jg 0 J5 ./Go 

rjn J 1 OO 

- C l JJ‘*i**«)\'*** + °l JjE a^ i) 

2 ^ r j n oq 

- G l J a \^*>)\^+cj s J G \E iXj{iX ) t+ 1 f ^-‘^ 


z 0J e iXjt ej 


dxdt 


dxdt 


(4.24) 


-5 ,=i '-5 


< C / / |iz(x, i; z 0 )| 2 cte<ii + (7 

J<5 Jg 0 

< C f I \tz(x,t; z 0 )\ 2 dxdt + C\A~^z 0 \y,. 

Js Jg 0 


'y ^2 j e,- 


dxdt 
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Step 2. In this step, we get rid of the term \A^Z{)\yi in the right hand side of (4.24) by 
a compactness-uniqueness argument. More precisely, we are going to prove that 


Ny' < C 


| tz(x, t; z 0 )\ 2 dxdt. 


(4.25) 


IS JGq 


If (4.25) is not true, then we can find a sequence {zq}™ = 1 C L 2 (G) with \zq \h- 2 (g) — 1 such 
that 

1 


/ / \tz(x, t; Zn)\ 2 dxdt < —. 

Is JGo n 


(4.26) 


Since {zq }™ =1 i s bounded in V', we can find a subsequence c { z o}ri=i suc h that 

ZQ k converges weakly to some z£ E H~ 2 (G) (4.27) 

and 


A^Zq* converges strongly to A^Zq in V'. 
According to (4.26) and (4.24), we know that 


(4.28) 


-Ink |2 


I^A 3 




0 ly ' - C n k 

This, together with (4.28), implies that there is a positive constant C > 0 such that 

1 


I^A z o\v' — 


C' 


(4.29) 


Thus, the limit Zq is non trivial. 

Further, from (4.21) and (4.22), we know that z(-, •; z^ k ) converges weakly to z(-, •; Zq) in 
L 2 (5,T-L 2 (G)). Hence, 


JS Jg 0 

Therefore, we find that 


z(x, t; ZQ)\ 2 dxdt < lim 

k—yoo 


1 S JGo 


\z(x, t; ZQ k )\ 2 dxdt < lim -— = 0. 

k^o o on k 


z(-, -;zq) = 0 in G 0 x (<5,T). 


(4.30) 


We would like to show that this leads to = 0 which would then yield to a contradiction. 
To do this, we introduce the linear subspace 


£ = {^o £ y' '■ the solution to (4.10) with the initial datum z 0 fulfills 
z(-, •; z 0 ) = 0 in G 0 x (S,T)}. 

Clearly, z^ given in (4.27) belongs to £. We want to prove that £ = {0}, which would be in 
contradiction with the fact that Zq is nonzero. 

Step 3. To show the claim that £ = {0} and conclude the proof, we proceed in several 
steps. 

Step 3.1. We first prove that £ C L 2 (G). 
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To do this, given £ e (0, 5) and any solution z, we introduce the discrete time-derivative 
4 (l , i;Zo) = !hlt±£i£h-£(£^ for te[0,T-S\. (4.31) 

Then we have that 
z £ (x,t]z 0 ) 


1 00 1 
ly-_-_ 

£ =1 ^A j(t + £ ) + 1 


z 0 ,je 


h+d e . 


1 A 1 


7 E 7 


£ i\d + 1 


z o,j e j e i 


3 = 1 


1 00 1 
±Y-. _-_ 

£ iXj (t+£) +1 


z 0J (e- a ^-e- a - t )e 1 + -J2 hr 


J 

OO 


£7 


z o,j- 




T' 1 i\ j (t +£)+1 UJ £ 


e -* A £e. 


-£77 


s . ^ LiAj (t-j-s) “I - 1 iAjt - ! - !- 

i\j 


z o,j e lXjt ej 


[i\j(t+£) + l](i\jt + 1) 


*oj-e iXit ej 


and 


Let 


z £ (x 


■ z ») = £ - £ 7 


j=i 


i=i 


i\j£ T 1 


z o,j£j- 


(4.32) 

(4.33) 


VAX,t) = ^X-7—, 


: z o,j- 


e -i\j£ _ 1 


e * A £ef 


^zAj(i + £) £ ' j~^iXj(t+£) 

Then, again, by assumption (A3), as in the proof of (4.23), we have 


£ 7T-7L7oV-“ j ‘<V (4.34) 


—4~ — — —£!: V =\z £ (x,0;zo)\ 2 v ,<C f I \{t + £)v E (x,t)\'dxdt, 
£ v ' Js Jg 0 

where the constant C is independent of £. Thus, we obtain that 


z £ (x,£;z 0 ) - z £ (x,0;z o ) 


£ 


V' 


pT—S p pT—5 p 

< C / / \z e {x, t; z 0 )\ 2 dxdt + C / / \z e {x, t; z 0 ) — v £ (x, t)\ 2 dxdt 

Js Jg 0 Js Jg 0 

pT—S p pT—S p 

< C / / \z £ (x, t; z 0 )\ 2 dxdt + C / / \z £ (x,t; z 0 ) — v £ (x,t)\ 2 dxdt. 

Js Jg 0 Js Jg 


(4.35) 


Let us estimate the second term in the right hand side of (4.35). From (4.32) and (4.34), we 
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have that 


cT—S 


15 jg 

pT—5 r , 00 


< 2 


1 5 JG 


z e (x,t;zo) — v £ (x,t) 'dxdt 

1 1 


^ Lx, 


3 = 1 

r T-S r , 00 


+2 


E 


j(t + e) + 1 iX j(t + e) 

iXj 


z o,j- 


e -i\j£ _ 1 


-e lXjt e, 


dxdt 


s J g ' [& A j (t e) T 1] (iXjt + 1) iXj (t + e) 


z o ,j e 3 e j 


dxdt 


If A j < 7 , then 


e -i\j£ _ ^ 


+ C J2j2 Zl 


\2 °,j’ 
j = 1 ^ 


e 


7lEiiHV)‘ = <vEsS2a,. 


If A,- > A, then 


£ I ^ /c! 

fc=l 


e -i\j£ _ ^ 


fc=l 


fc=l 


k\ 


= -\e 
£ 


—i\,£ 


l| < Xj\e~ ix i £ - l| < 2Xj. 


According to (4.36), (4.37) and (4.38), we obtain that 

n 

J 5 Jg 


z e (x,t ; z 0 ) - v £ (x,t)\ 2 dxdt < j2 z o,j- 

j =1 J 


As a result of (4.35) and (4.39), we have that 


z £ (x,E]z 0 ) - z £ (x,0;z 0 ) 


< C 


V' 


cT-5 


is Jg 0 


|A £ (x,t; z 0 )| dxdt + T^^o.r 

3=1 i 


Since z 0 E £, we know that 

4(x, t; z 0 ) = 0 in G 0 x (d, T). 
This, together with (4.40), implies that for any e G (0,d), 


z £ (x,£;z 0 ) - z £ (x,0;z o ) 


00 9 

zi 


V’ 


<cEA<c np 

j-i 


Letting £ —>■ 0 in (4.41), we obtain that 


00 ~2 

•^n 


ko|ia (G) < C \J2 iX i z °J e i v , - - C'l^olv'- 

i=i 1=1 i 


Then, we have that S C L 2 (G). 
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(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 

(4.42) 






























Step 3.2. The same estimate deduces that £ is a finite-dimensional subspace of L 2 (G), 
since 

\ z o\l 2 (G) — C\zq\y,. (4.43) 

Step 3.3. We now claim that A&£ C £. 

Utilizing (4.21) again, we have that 

z(-, •; zq) E C([0, T];L 2 (G)) for any z 0 £ £. 

Since z(-, •; z 0 ) = 0 in G 0 x (S, T ) for any S > 0, we see that z 0 = 0 in G 0 . 

Thanks to £ C L 2 (G), we get that A&£ C V' and z(x, t; A&zo) = A&z{x,t\ zq) = 0 in 
Go x (5, T). Therefore, we obtain that A^£ C £. 

Step 3.4. To conclude, assume that £ ^ {0}. Then, there would exist a non-trivial 
eigenfunction ip £ £ and an eigenvalue /i£l such that 

{ —= flip in G, 

ip = 0 on dG, 

ip = 0 in Go- 

However, by the classical unique continuation property for elliptic equations, ip = 0 in G, 
which contradicts that ip is an eigenfunction. 

This concludes the proof of the fact that £ = {0}. 

As a consequence, we derive the desired estimate (4.20). □ 

Remark 4.7 We have utilized a compactness-uniqueness argument in the above proof, which 
has been extensively used in the proof of observability estimates (see [7] for example). Note 
however that, normally, this is done for solutions of PDE models. The averages under 
consideration not being solutions of a specific PDE this argument needs to be carefully adapted 
as we have done above. 


5 Further comments and open problems 


5.1 Further comments 

5.1.1 Boundary averaged control for random heat equations 


In this subsection, for convenience, we assume that dG is C°° smooth, although most 
comments and results make sense with different weaker regularity assumptions. 

We have solved the internal averaged controllability problems for some particular classes 
of random heat and Schrodinger equations. The same could be done for boundary control 
problems. 

Let us consider the following heat equation with boundary control and random constant 
diffusivity: 


Vt ~ o:Ay = 0 
y = u 

y = 0 

, 1 /( 0 ) = Vo 


in G x (0, T\, 
on T 0 x (0,T), 
on (<9G\T 0 ) x (0,T), 
in G. 
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(5.1) 




Here r 0 is an open subset of dG, a(-) is a random variable, u G L 2 (0, T; L 2 (T 0 )) and yo G 
L\G). 

We have the following result. 


Theorem 5.1 Let a(-) be a uniformly distributed or an exponentially distributed random 
variable. The system (5.1) is null controllable in average with control u(-) G L 2 (0, T; L 2 (T 0 )). 
Further, there is a constant C > 0 such that 


Ml 2 (0,T;L 2 (To)) < C|2/o|l 2 (G)- 

To prove Theorem 5.1, we consider the adjoint system of (5.1) as follows: 

= 0 in G x (0,T], 
z — 0 on dG x (0, T), 

z(T) = zo in G. 

One only need to prove the following result. 


(5.2) 


(5.3) 


Theorem 5.2 There exists a constant C > 0 such that for any zo G L 2 {G), and either ck(-) 
is a uniformly distributed or an exponentially distributed random variable, it holds that 


z(-, 0,(n; z 0 )dF(u) 


L 2 (G) 


< c 


' o Jr o 


dz(x, t, u>] zq) 
dv 


dP(o;) 


dTodt. 


(5.4) 


The proof is very similar to the one for Theorem 3.2. We only give a sketch here. Let us 
assume that a(-) is an exponentially distributed random variable. From [34, Page 345], we 
have the following result: 


^2 a 2 < C 2 e c *^. 


c t 2~ t l 


f / 

hi J r 0 


de^ 2 


V evh f a ._i dTodt for any 0 < t x < t 2 < T. (5.5) 
dv 


Let T fc = (1 — )T for k G N and r k = 2 2 ( fc+1 )[ln(6C2) + C 2 \. Similar to the proof of 
(3.25), we can obtain that 


Tk+i~T k \ z( rr . „ w C 2 e C 2 ^+l {rr ^ v 

QC 2 e C2 ^ Z ^' ,Ik,A ' 0L2< ' G) C 2 e C2 '/ ¥ * ifc+1 Ikt 


L k+l~ 1 k 


rk 


\z(-,T k+1 -, ^o)|l 2 (g) 


'fc+i 


< 


T 0 


dz(x,t] z 0 ) 


dv 


dT 0 dt. 


< 


By summarizing the inequality (5.6) from k — 1 to k — oo, we obtain that 

f T r dz(x,t-,zo) 

k =i ^ r ° 

where 

T k +2 — T k + 1 C 2 e C2 '^ 7 + 1 m v T k+i~ i k 


rj~\ rjn TXD 

6(72e C2 v^^°^ 2(G) + S Tfc+i; z ° )I L2 (G) 


dv 


dTodt , 


(5.6) 


(5.7) 


fk = 


6G 2 e C2 ^Tl C 2 e C2 'TT 


■(Tk +i — T k )e 
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rk 


, k = 1,2,-•• 


2 
























From (3.18) and r k = 2 2 ( fc+1 ) [ 111 ( 662 ) + C* 2 ], we have that 


fk > 0 for any k = 1 , 2 , 
This, together with (5.7), deduces that 


\ Z o\L 2 (G) - 

This completes the proof. 


6 C 2 e c ^ r T r 

J 0 J To 


t 2 -t 1 


dz(x,t;z 0 ) 


dv 


dT 0 dt. 


(5.8) 

□ 


5.1.2 Boundary averaged control for random Schrodinger equations 

Consider the following random Schrodinger equations: 

yt — aiAy = 0 in G x (0, T], 
y = u 

y = 0 

, 2 /( 0 ) = Vo 

Here T 0 is an open subset of dG, a(-) is a random variable, u G L 2 (0, T; L 2 (T 0 )) and y 0 G 
L 2 {G). We have the following controllability results. 

Theorem 5.3 System (5.9) is null controllable in average if a is a random variable with 
normal distribution or Cauchy distribution. 


on T 0 x (0,T), 
on (dG\T 0 ) x ( 0 ,T), 

in G. 


Further, we assume the following condition holds: 

(A4) Whatever T > 0 and k 7 ^ 0 are, there are constants C 3 and C 4 such that for any 
tpo G Hq(G), the solution </?(-, •) to (4.13) satisfies 


W 


Olifi(G) ^ C 3 


'Go 


dip 


dv 


2 

dYftdt ft C41 I (g?) * 


(5.10) 


We refer the readers to [33] for the conditions on Tq for which (A4) holds. 


Theorem 5.4 Assume that (A4) holds. We have the following results: 

• Let «(•) be a uniformly distributed random variable on an interval [a,b\. Then, the 
system (4.1) is exactly controllable in average with V = H = L 2 (G ) and U = i7 _1 (r 0 ). 

• Let «(•) be an exponentially distributed random variable. Then, the system (4.1) is 
exactly controllable in average with H = L 2 (G), V = Hq(G) and U = L 2 (r 0 ). 

• Let «(•) be a random variable with Laplace distribution. Then, the system (4.1) is 
exactly controllable in average with H = L 2 (G), V = H 3 (G ) HH^G) and U = L 2 (T 0 ). 

• Let «(•) be a random variable with standard Chi-squared distribution. Then, the system 
(4.1) is exactly controllable in average with H = L 2 (G), V = H k ~ l {G) fl Hq{G) and 

U — I/ 2 (r 0 ). 











As usual, we introduce the adjoint system of (5.9) as follows: 


{ z t + iaAz = 0 in G x (0, T), 

z = 0 on 9G x (0,T), (5-11) 

z(T) = zq in G. 

We can prove the following results. 

Theorem 5.5 Let a(-) be a standard normally distributed random variable. There exists a 
constant C > 0 such that for any zo G L 2 (G), it holds that 


E 

3 = 1 


Znje 


-A fT 2 


L 2 {G) 


< c 


To 


dz(x, t, oj] zo) 
dv 


dP(u) 


dTodt. 


(5.12) 


Theorem 5.6 Let a(-) be a random variable with standard Cauchy distribution. There exists 
a positive constant C such that for any z 0 G L 2 (G), it holds that 


J2 z °’i e XjTe i 

3 = 1 



dz(x, t, u >; zo) 
dv 


dP(u) 


dTodt. 


(5.13) 


By means of (4.12), Theorem 5.6 is nothing but the boundary observability estimate for 
heat equations. The proof of Theorem 5.5 is very similar to the one for Theorem 5.2. We 
omit it here. 

Further, the proof of Theorem 5.4 is also very analogous to the proofs for Theorems 4.3. 
We only give a sketch of the proof of the second conclusion in Theorem 5.4. 

Let z 0 G Lf" 1 (G) in (4.10) and z(x, f; £ 0 ) = f Q z(x, T — t, cu; z 0 )dF(uj). Then, we only need 
to prove that 


z 0 


H-HG) 



dz(x,t,z 0 ) 

dv 


2 

dT 0 dt. 


(5.14) 


Assume that zq = z o,j e j- We have 


z(x,t] z 0 ) 


OO 


E 


1 

iXjt + 



(5.15) 


Let 

v(x,t) = J2fx^ z °’^ iXjte r 

3 = 1 j 

From (A4), for a fixed 6 G (0, T), we have that 


zo 


H-^G) 


E 


3 = 1 



e * Ajt e 


2 

HfcG) 



t 


dv(x, t) 
dv 


2 

dT 0 dt. 


(5.16) 
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Therefore, 


Folff-i^G) 

pT p 


< c 

< c 

< c 

< c 


>s J r 0 


dv(x,t) 


dv 


dT 0 dt 


\f [ 

dz(x,t-z 0 ) 

2 dr 0 dt + [ [ 

+ dv(x,t) + dz(x,t-,z 0 ) 

2 1 

dT 0 dt 

L Js J r 0 

dv 

Js J r 0 

dv dv 



(5.17) 


Is J r 0 

pT p 


dz(x,t;z 0 ) 


r T oo 


dv 

dz(x,t-,z 0 ) 2 


dT 0 dt + C 


' 7 [\ (i \ -1 


Js J r 0 

We claim that 


<9i/ 


Is I iXjiiXjt + 1 ) 


■soj-e ^ej 


HHG) 


dt 


\ z o\‘h- 1 (G) — C 


' 0 J Tq 


<95(x,t;z 0 ) 


<9z/ 


dT 0 dt. 


(5.18) 


If (5.18) is not true, then we can find a sequence {zq}™ =1 C L 2 (G) with |zq \h~ 2 (g) — 1 such 
that 

r T r t dz(x,t-,zo) 2 . 1 


' 0 V To 


<9^ 


dTodt < —. 
n 


(5.19) 


Since {^g}^=i is bounded in H 1 (G), we can find a subsequence c {^ 0)^=1 such 

that ZQ k converges weakly to some Zq G H~ l {G). From (5.19), we know that 


12 


1 

> — 


Z 0\H~3(G) — £? 


for a positive constant and 


dz(-,-; z o) 

dv 


= 0 on r 0 x ( 0 , T ). 


(5.20) 


(5.21) 


Put 


£ = {x;o G H 1 (G) : the solution to (5.11) with the final datum Zo fulfills (5.21)}. 

Analogous to the proof that £ = {0}, we can prove that £ = {0}, which contradicts (5.20). 
Hence, we know that (5.18) holds. □ 

5.1.3 Averaged control for random heat and Schrodinger equations from mea¬ 
surable sets 

We have considered the averaged control problems of random heat and Schrodinger equa¬ 
tions for the internal control (resp. boundary control) supported in Go x E (resp. To x E), 
where Go C G (resp. r 0 C dG) is a nonempty open subset. By means of the method in [4], 
one can consider the case that the internal control (resp. the boundary control) is supported 
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in a measurable subset PcGx(0,T) ( resp. Q C dG x (0, T)). For instance, we can consider 
the following systems: 


and 


{ Vt - aA y = xvu 

y = o 

2/(0) = 2/o 


in G x (0, T), 
on dG x (0, T), 
in G, 


{ Vt - iaAy = xvu 
2 / = 0 
2/(0) = 2/o 


in G x (0, T), 
on dG x (0, T), 
in G. 


(5.22) 


(5.23) 


Here y 0 G L 2 (G), X? C G x (0, T) is a Lebesgue measurable set with positive Lebesgue 
measure and u G L°°(G x (0, T)). One can combine the proof of Theorem 3.1 and Corollary 
1 in [4] to prove the following results. 


Theorem 5.7 Assume that, either a(-) is a uniformly or exponentially distributed random 
variable. Then, the system (5.22) is null controllable in average with control u(-) G L°°(G x 
(0, T)). Further, there is a constant C > 0 such that 

Ml°°(Gx( 0,T)) < C|2/oU 2 (G)- 


Theorem 5.8 If a is a random variable with normal distribution or Cauchy distribution, 
then the system (5.23) is null controllable in average with control u(-) G L°°(G x (0,T)). 
Further, there is a constant C > 0 such that 

Ml°°(Gx(o,t)) < (■ !2/o|l 2 (G)- 

Next, let us consider the following systems: 

2 h ~ ®Ay = 0 
y = u 

y = o 

, 2 /( 0 ) = 2/0 

and 

y t — aiAy = 0 
y = u 

y = o 

, 2 /( 0 ) = 2/0 

Here Q is a Lebesgue measurable subset of dG x (0,T) with positive Lebesgue measure, «(•) 
is a random variable, u G L°°(dG x (0, T)) and y 0 G L 2 (G). 

One can combine the proof of Theorem 5.1 and Corollary 1 in [4] to prove the following 
results. 


in G x (0, T ], 
on Q, 

on [dGx (0,T)]\g, 
in G, 


(5.24) 


in G x (0, T], 
on Q, 

on [dG x (0, T)] \ Q, 

in G. 


(5.25) 
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Theorem 5.9 Assume that, either a(-) is a uniformly or exponentially distributed random 
variable. Then, the system (5.24) is null controllable in average with control u(-) G L°°(G x 
(0, T)). Further, there is a constant C > 0 such that 

Ml°°(9Gx(0,T)) < C\yo\L 2 (G)- 

Theorem 5.10 If a is a random variable with normal distribution or Cauchy distribution, 
then the system (5.25) is null controllable in average with control u(-) G L°°(G x (0,T)). 
Further, there is a constant C > 0 such that 

Mi°°(dGx(o,T)) < C\y 0 \L 2 (G)- 


5.1.4 Internal averaged control for random fractional Schrodinger equations 

We have studied the averaged controllability problems for some random heat equations 
and random Schrodinger equations. In the results proved so far we have obtained averaged 
controllability for parameter-depending equations that were controllable for each value of 
the parameter. Here, we give an example of model which is not null controllable when one 
fixes an u but that gains null controllability by the averaging process. 

Consider the following equation: 

iiy t + aA\y = Bu in(0,T], 

1 2/(0) = 2/o- ( ' } 


Here 7 G (|, |), y 0 G L 2 (G), u G L 2 (0, T; L 2 (G 0 )) and Bu = y Go u. 
The adjoint system of (5.26) reads 


iz t — aA\z = 0 in [0, T), 
z{T) = Zq, 


(5.27) 


where z 0 G L 2 (G). Assume that z 0 = Y2'jLi z o,j e j- If Q; (') : H —> M is a standard normally 
distributed random variable, then we know that 


z(x, t, u>] z 0 )dF(co) 



OO 


= 

3= 1 


A^(T-t) 2 


-'3- 


(5.28) 


Similar to the proof of Theorem 4.1, we can establish the following result. 

Theorem 5.11 Let E C [0, T] be a measurable set with positive Lebesgue measure m(E). 
There exists a constant C > 0 such that for any z 0 G L 2 (G), it holds that 


E 

3 = 1 


Zoje 


-A CT 2 , 


L 2 (G) 


< c 


z(x, t, uj; z 0 )dF(uj) 


L 2 (G 0 ) 


dt. 


(5.29) 
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Theorem 5.11 implies that the system (5.26) is null controllable in average. However, it 
is not null controllable for any fixed a G R even for d — 1. For example, let a = 1 and 
G = (0,1). Then the solution to the system (5.27) reads 


z(x,t) = zp^e h-hd 27 * sm(jnx). 

3 = 1 


g mC6 

lim {[(j + l)vr)] 27 - (jtt) 27 } = 0, 

o 

we know that the following inequality 


3 = 1 


Z 0,3 < C 


| z(x, t)\ 2 dxdt 


>G 0 


does not hold for any C > 0. 


5.1.5 Internal averaged control for random heat equations with variable coeffi¬ 
cients 

We can also consider the approximate averaged controllability problem of some more 
general random heat equations. We make the following assumptions on the coefficients 
a jk : Gxll4 M nxn (j, k = 1, 2, ■ • • , n): 

(HI) a? k {-,u) : G —* M is analytic, W-a.s., and a^ k = a k F 

(H2) For a.e. ui G there is a constant C{uj) > 0 such that for any multi-index 
V= (vi,--- ,Vn) e (NU{0}) n , 


d v a jk (x,u}) C(uj)\r]\ 


dx^ 


< 


m 


for any j,k = !,■■• ,n, 


where R is a positive constant larger than max xg g \x\, and C{-) satisfies that 

/ C(oj)dF(u3) < oo. 


(H3) There exists a constant So > 0 such that 

n 

Y t, x)£ J £ k > so|C| 2 , v x, 0 = (ujjXjg 1 , ■ ■ • ,£ n ) efl x G x R n . (5.30) 

j,k=1 

Consider the following random heat equation: 

' n 

yt - Y ( a3k y*i) Xh = XG ° U in G x (°> T )’ 

< j ’ k=1 (5.31) 

y = 0 on dG x (0, T), 

k y{ o) = yo in G. 

Here the initial datum y 0 G L 2 (G). 

We have the following result. 
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Theorem 5.12 Under the assumptions (H1)-(H3) above system (5.31) is approximately 
controllable in average in any time T > 0 and from any open non-empty subset Gq of G. 


Proof: By Theorem A.3, we only need to prove that the adjoint system of (5.31) satisfies 
the averaged unique continuation property. Its adjoint system reads 

f n 

A + ( a 3k z X] ) Xk = 0 in G x (0, T), 

< j ’ k=1 (5 321 

z = 0 on dG x (0,T), 

z(T) = zq in G, 


where the final datum zo £ L 2 (G). From (HI) and (H3), we know that for any t £ [0, T) 
and a.e. cu £ fl, the solution z(-,t,uj;z o) is analytic in G (see [18, 28] for example). Further, 
for any ball B r C G with radius r, there is a constant C > 0 such that for any multi-index 
7 ] £ (NU{0}) n , 


d v z(-, t, Ui; Zq) 

dxv 


< CC(uj)- 


in B r . 


From (H2), we have that 

d v J n z(-,t,u)dP(u) <C W 
dx 71 ~ rid 

Then, we know that f n z(-, t, oj)dF{uj) is analytic in B r . Hence, it is analytic in G. Since 
f Q z(-, -,uj)dP(cj) = 0 in G 0 x (0, T), we get that it vanishes everywhere in G x (0, T ). Noting 
that it is continuous in Li 2 {G) with respect to t, we conclude that Zq = 0 in G, which implies 
that (5.32) satisfies the averaged unique continuation property. □ 


5.1.6 Averaged controllability problems for the random heat and the random 
Schrodinger equations random initial data 

One can also consider the internal and boundary averaged controllability problems for 
random heat and random Schrodinger equations with random initial data. Let us first 
consider the following random heat equation: 


{ y t - a Ay = Xg 0 xeu in G x (0, T), 

y = 0 on 9G x (0,T), (5.33) 

y(0,uj) = y 0 (uj) in G, 


Here y 0 (•) £ L 2 (12; L 2 (G)), and G 0 and E are subsets of G and [0, T] respectively, where the 
controls are being applied. The constant diffusivity a : 12 —?• M + is assumed to be a random 
variable. 

According to Remark A.l, we know that to prove the averaged null controllability of 
(5.33), we only need to establish the following observability estimate: 


x(0,o;;zo) 


dP(u) ) <C 

LUG) V J E 


z(t, c u; z 0 )^IP( w ) 


LHGo) 


dt , 


(5.34) 


where z solves (3.7) and C is independent of 2 0 . 
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If «(•) is the exponentially distributed random variable, then 



£(0,o;;zo) 


d¥(oj) 

L 2 (G) 



E 

3 =1 


_I_~2 —2\jT 

2XjT + 1 0j 


(5.35) 


From (5.35), similar to the proof of the inequality (3.8), one can obtain (5.34). The same 
thing can be done if a(-) is a uniformly distributed random variable on [a, b] for 0 < a < b. 


Proposition 5.1 Let (Al) and (A2) hold. Assume that, either a(-) is a uniformly or 
exponentially distributed random variable. Then, the system (5.33) is null controllable in 
average with control u(-) G L 2 (0, T; L 2 (Go)). Further, there is a constant C > 0 such that 

Ml 2 (0,T;L 2 (Go)) < C|2/o|l 2 (0;L 2 (G))- (5.36) 

Thanks to Remark A.l, we know that in order to show that (5.33) is approximately 
controllable in average, one just needs to prove that the solution to (3.7) satisfies the averaged 
unique continuation property, which is obtained in the proof of Theorem 3.3. Hence, we have 
the following result: 

Proposition 5.2 Let (Al) and (A2) hold. System (5.33) is approximately controllable in 
average, provided that cc(-) is a uniformly distributed or an exponentially distributed random 
variable. 


Next, we consider the averaged controllability problem for the following random 
Schrodinger equation: 

{ y t - iaAy = xg 0 xeu in G x (0, T), 

y = 0 on dG x (0, T), (5.37) 

1/(0) = y 0 in G. 

Here yo(-) G L 2 (Ll]V), and Go and E are subsets of G and [0, T] respectively, where the 
controls are being applied, a : fl —> M is assumed to be a random variable. 

In virtue of Remark A.l, we know that to get the averaged exact controllability of the 
system (5.37), one only need to prove that the solution to the equation (4.10) is exactly 
averaged observable. As a result of these facts, we know that the conclusions in Theorem 
4.3 also hold for the system (5.37). More precisely, we have the following results: 

Proposition 5.3 The following results hold: 

• Let «(•) be a uniformly distributed random variable on an interval [a, b\. Then, the sys¬ 
tem (5.37) is exactly controllable in average with V = H = L 2 (G) and U = H~ 2 (G 0 ). 

• Let a(-) be an exponentially distributed random variable. Then, the system (5.37) is 
exactly controllable in average with H = L 2 (G), V = H 2 (G) fl H( (G) and U = L 2 (Gq). 
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• Let a(-) be a random variable with Laplace distribution. Then, the system (5.37) is 
exactly controllable in average with H = L 2 (G), V = H 4 (G) D Hq (G) and U = L 2 (G 0 ). 

• Let a(-) be a random variable with standard Chi-squared distribution. Then, the system 
(5.37) is exactly controllable in average with H = L 2 (G), V = H k (G ) D Hq(G) and 
U = L 2 (G 0 ). 

Further, let us consider the averaged null controllability problem for the system (5.37). 
Due to Remark A.l, we only need to prove the following observability estimate: 


z(Q,lv,zq) dPhu)) < C / 

l 2 (G) V J E 


z(x, t, l u; £o)dIP(a;) 


L 2 (G 0 ) 


dt. 


When a(-) is a normally distributed random variable, we have that 


z(0, u;;2 0 ) dP(u;) = —__ , 

LHG) y/2^ 


/ OO 2 °° 

e~^Y, z li e ~ 2iXjaTda 

•°° j=i 




2 -2A ?t 2 
0 


j =1 


(5.38) 


(5.39) 


When a(-) is a random variable with the Cauchy distribution, we have that 


z(0, ui; Zq) d¥(co) = 
L 2 (G) 


.Loo ?r(l + a 2 ) 

OO 


2iXjaT da 


■e~ 2XjT . 


o =i 


(5.40) 


Thanks to (5.39) and (5.39), and similar to the proof of Theorem 3.2, one can show 
that the inequality (5.38) holds when a(-) is a random variable with normal or Cauchy 
distribution. Therefore, we have the following result: 


Proposition 5.4 Under the assumptions of Theorem f.l the system (5.37) is null control¬ 
lable in average if a is a random variable with normal or Cauchy distribution. 


5.1.7 Control of the variance of the state of the system (2.1) 

If a system is null( resp. exactly) controllable in average, then we can drive the average of 
the state to rest (resp. a given destination). It is then natural to analyse whether one can 
control the higher order moments of the state and, in particular, the covariance Cor(y) of 
the state y, which is defined as follows: 

Cor (y) = I (y(T, uj) - J y(T,iv)d¥(uj )) ® (y{T,w) - j y(T,uj)dF(oj)^jdF(uj), 
where <g) denotes the tensor product of two elements in H. 
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Cor(y) is an element in H® H, which is also a Hilbert space. With the goal of reinforcing 
the averaged controllability property studied in this paper, a first thought could be to look 
for the controllability of the covariance. But this is not suitable since the variance measures 
how far the random variable is spread out. In other words, loosely speaking, the smaller 
the \Cor(y(T))\ H(S)H is, the closer the state y(T, •) is concentrated to an element in H and, 
therefore, closer to the simultaneous controllability property. A more natural goal is to look 
for the null controllability of the covariance. However, from the definition of Cor(-), it is 
trivial to see that Coi(y(T)) = 0 if and only if y(T , •) = J Q y{T,uj)dF(uj) with probability 
one, i.e., y(T,u>) is independent of u> with probability one. Nevertheless, this is impossible 
if we only use parameter independent control (see Remark 2.4). Therefore, since we want 
to make \Cor(y(T))\v®v as small as possible, it is natural to study the following optimal 
control problem: 

Problem (OP2): Minimize |Cor(|/(T))|y ( g ) y for 

u(-) G U nu u[0,T} = j?i G L 2 (0,T; U ) : the average of solutions to (2.1) corresponds 
to u satisfies that 

Problem (OP2) is an optimal control problem with a terminal constraint. Following [36, 53], 
one can derive a Pontryagin type Maximum Principle for it. However, this is beyond the 
scope of this paper. 

5.2 Open problems 

There are many interesting and important (at least for us) problems in this topic. We present 
some of them here briefly: 

• Averaged controllability problems for general random heat and Schrodinger 
equations. 

We have only studied the averaged controllability problems for some very special classes 
of random heat and Schrodinger equations, for which the averaged dynamics could be 
computed explicitly. It would be interesting to investigate some more general classes 
of parameter dependent systems. For example, is it (3.1) null controllable in average 
if we take a to be a random variable with Chi-squared distribution? 

Furthermore, the method we use to study (3.1) and (4.1) depends on the fact that the 
eigenfunctions of — aA/\ are independent of ui. Thus, more general random heat and 
Schrodinger equations (as, for instance, (5.31)) where the principal part of the PDE 
depends on co cannot be treated in this way. 

The method of proof employed to show the approximate averaged controllability of 
(5.31) is based on the use of the space-time analyticity properties of solutions (5.32) 
to derive unique continuation properties, but it does not provide any quantitative 
information. 


y(T,cj;y 0 )du = 0 
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The relationship between averaged controllability properties and the ran¬ 
dom variable. 

We have shown that different random variable a(-) may lead to different controllability 
property of the system (4.1). It is interesting to give a description of the relationship 
between the random variable a(-) and the controllability property of the system (4.1). 
For instance, for what kind of random variables, the system (4.1) is null controllable 
in average? Is there a random variable such that the system (4.1) is neither exactly 
controllable in average nor null controllable in average? 

The null averaged controllability problem for random fractional heat equa¬ 
tions. 

We have proven that the random fractional Schrodinger equations are null control¬ 
lable in average for 7 £ (|, 1). It is more interesting to study the same problem for 
random fractional heat equations, which describe the anomalous diffusion process. In 
particularly, consider the following system: 


y t + aA\y = Bu in( 0 ,T], 
1/(0) = Vo- 


(5.41) 


Here 7 > 0, yo E L 2 (G), u E L 2 (0, T; L 2 (Gq)) and Bu = Xg 0 u - 

It is clear that when 7 < 1 and a(-) is a uniformly or exponentially distributed random 
variables, the system (5.41) is not null controllable in average. However, is it possible 
to find some random variable a(-) such that the system (5.41) is null controllable in 
average for some 7 e (0,1]? 

Averaged controllability problems for general random evolution partial dif¬ 
ferential equations. 

The method used in this paper can only be applied to the case of equations in which 
the random operator is aA, which is very restrictive. 

Several methods have been developed to solve the controllability problems for deter¬ 
ministic partial differential equations. For the heat equation, the existing methods 
include Carleman estimates ([23]) and the moment method ([5]). For Schrodinger 
equations, Carleman estimates can also be applied ([30]), together with the moment 
method ([5]), multipliers ([42]), microlocal analysis ([32]), etc. 

It would be interesting to generalize these method to deal with the averaged control¬ 
lability problems for more general random heat and Schrodinger equations. The paper 
[2] contains an interesting survey on the controllability of parabolic systems. It would 
be interesting to explore possible applications to averaged control. 


Averaged controllability problems for nonlinear random evolution equa¬ 
tions. 
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We have studied the averaged controllability problems for linear random evolution 
equations. The same problem can be considered for nonlinear random evolution equa¬ 
tions. A possible method to handle the nonlinear problem is to follow what people do 
for the classical controllability problems, that is, combining the controllability result 
and some fixed point theorem or inverse mapping theorem. However, as the average 
of the solution of linear transport equation with respect to velocity helps people study 
the nonlinear transport equations, we expect that one can get better result than the 
ones obtained by applying the method mentioned above directly. For example, can 
one prove that a random Schrodinger equation with a cubic term is exactly or null 
controllable in average? 

• Numerical approximation of averaged controls. 

The numerical approximation of control problems is studied extensively in the litera¬ 
ture. We refer the readers to [15, 20, 55] and the rich references therein for this topic. 
This question also arises in the context of averaged control. This can be done based 
on the variational characterization of the average controls given in this paper. 


A Appendix 


A.l Reduction of averaged controllability to averaged observabil¬ 
ity 

We have the following results. 

Theorem A.l System (2.1) is exactly controllable in average in E with the control cost 
C > 0 if and only if the adjoint system (2.8) is exactly observable in average in E. 

Proof of Theorem A.l: The “if 7 part. Let us define a linear subspace X C L 2 (E] U ) as 

X = jxe(-) J B*(uj)z(-,uj-Zo)dF(uj) : z 0 E V' 

and a linear functional A on A as 

f(xe{-) J B*(uj)z(-,u-,z 0 )d¥(uj )) = {yi, z 0 ) V y - (y 0 , / z(0,u; z 0 )d¥(uj) 

From (2.9), we know that 

f(xe(-) f B*(u)z(-,u]Zo)¥(u) 

J 

< C(\y 0 \ v + \yi\v) (\zq\v + [ z(0,u}] z 0 )dP(cj) 


vy 


< C(\y 0 \ v + \yi\v) Xe(-) / B*(u)z(-,u; zq)¥(u) 


V' 


L 2 (0,T\U) 


Hence, F is a bounded linear functional on X with norm < C(\y 0 \v + \yi\v)- Then, 

it can be extended to a bounded linear functional on E 2 (E; U) with the same norm. We still 
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denote by F the extension if there is no confusion. Then, by Riesz Representation Theorem, 
there is a u(-) G L 2 (E ; U) such that for any v(-) G L 2 (E ; U), 

F(v) = (v,u) L 2 ( E . v) 


and 


u\ L 2 (E-,U) — |-F|£(*,]R) < C(\y 0 \v + \yi\v)- 


From the definition of F(-), we know that for any Zq G V', 


(yi,z 0 ) V y- (yo, [z(0,(j]Zo)dP(u)\ = [ XE(t)(u(t), I B*(u)z(t,u-,Zo)dP{u)\ dt. 

' Jn 1 vy Jo N Jn ' u 

(A.l) 

We claim that u(-) is the control we need. Indeed, taking the dual product of V' and V of 
z = z(t,u ; Zq) with ( 2 . 1 ) and integrating the product with respect to t in ( 0 , T ) and c o in 0 , 
we obtain that 


/ XE(t)(u(t), / B*(u)z(t,u;z 0 )dF(u)) dt 

= [ Xs(t) [ (u(t), B*(u)z(t,u] z 0 ))udF(uj)dt 
Jo Jn 

= [ Xe(J) [ (B(u)u(t), z(t, uj; zo)) v ,v'dP(uj)dt 
Jo Jn 

= / (y(T,u};y 0 ),zo)v,v'dF(v) - / (yo, z(0, cu; z 0 ))vydF(uj) 


/ y{T,u-,y Q )<JF(u),zo) - (yo, / z(Q,w,zo)dP(u) 

In /vy \ Jn ivy 


From (A.l) and (A.2), we conclude that 


y(T,uj]y 0 )dF(uj),z 0 \ , = {y 1 ,z 0 )vy, V ^ 0 e V', 


which deduces that f Q y(T , oj\ y 0 )dF(u>) = y x . 

The “only if’ part. Let zq G V'. We choose yo,yi G V' which satisfy that 

\yo\v = \yi\v < 2, 

/ z( 0 , a;; 2 0 )dP(a;) 


— ( 2 / 0 , / £( 0 ,w;z 0 )dP(w) 

\ Jn / W 

(yi, zo)vy = |In'- 


'O 


V' 


Let w(-) be the control such that 

M L 2 {E;U) < C(|yo|y + \yi\v) < C 


(A. 2 ) 


(A.3) 


and 


I y{T,u>iy 0 )dP(u>) = yi. 
n 
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Then, from (A.2), we have that 


\zo\v' + 


z(0,u;; zo)dF(cu) 


= / XE(t)(u(t), / B*(u)z(t,u;zo)dF(u})) dt. 

v ' Jo x Jn ' u 


Thus, we find that 


po \v 


2(0,0;; 2 0 )dIP(o;) 


V' 


< ML 2 (£;t/) 


_B*(o;) 2(-, cu; 2 o)gOP(u;) 


L 2 {E\U) 


This, together with (A.3), implies that 


|^o|v' — C 


B*(co)z(-,ur, 2 0 )oflP(u;) 


L 2 (E;U) 


Theorem A.2 System (2.1) is null controllable in average in E with control the cost C > 0 
if and only the adjoint system (2.8) is null observable in average. 

The proof of Theorem A.2 is very similar to the one for Theorem A.l. We omit it here. 

Theorem A.3 System (2.1) is approximately controllable in average in E if and only if the 
adjoint system (2.8) satisfies the averaged unique continuation property in E. 

Proof of Theorem A.3: The “if’ part. Since the system (2.1) is linear, we may assume that 
r/o = 0. Then, we only need to prove the following set 


At = 


y(T, iv; O)dP(ca) : y solves (2.1) with some control u( 


is dense in V. We do this by contradiction argument. If At was not dense in V, then we 
can find a tp £ V' with \<p\v = 1 such that 

(-0, (p)vy =0, V0 e At- 

On the other hand, similar to (A.2), we have that 


/ XE(t)(u(t), / B*(uj)z(t,uj-,z 0 )dF(u )) dt=( / y(T,u]0)dF(u),z o 

Jo x Jo. ' u x Jn ' v y 

Let z 0 = in (A.4). We have that 

f XE(t)(u(t), ( B*(uj)z(t,uj-,(p)dF(uj)^ dt = 0, Vw(-) G L 2 (0, T; U). 

Hence, we find that 


(A.4) 


Xe(-) / B*(u)z(-,u;<p)dP(u) = 0 in L 2 (0, T; U), 
Jn 

which implies that ip = 0 and leads to a contradiction. 
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The “only if’ part. We utilize the contradiction argument again. We assume that there 
is a z 0 E V' with \zq\v> = 1, such that 

Xe(-) [ B*(uj)z(-,uj;z 0 )dP(uj) =0 in L 2 (0, T; U). 

Jn 

This, together with (A.4), implies that 


( [ y(T,uj;0)dP{uj),z o ) =0, Vn(-) E L\E; U). 

' Jo 1 vy 


(A.5) 


On the other hand, since (2.1) is approximately controllable in average, we can find a u(-) E 
L 2 (E; U ) such that 

[ y(T,uj;0)d¥(uj) - z 0 


It is clear that for this f Q y(T, uj; 0)dP(u;), 

f y(T,uj;0)dP(uj),z o 


vy < 2’ 


1 

> 

vy 2 


which contradicts (A.5). 


□ 


Remark A.l One can also consider the case where yo G L 2 (f2; V), i.e. when the datum to 
be controlled depends on the parameter oj as well. 

The proof of Theorem A.l applies replacing the terms (y 0 , f Q z(0,u; z 0 )dP(u;))vy and 
\Vo\v by the terms f Q (y 0 (ca), z(0, uj; Zo))vydF(uj) and \yo(-)\L 2 (n:V)> respectively. The same 
can be said about Theorem A.3. 

The situation is different and much more delicate in the context of averaged null con¬ 
trollability. Note that, when considering initial data to be controlled depending on uj, one 
requires an observability estimate of the form 


z(0,ui;zo) d¥(u) < C / XE(t) / B(ui)*z(t, u>; zo)dF(co) 

V' .In .lo 


dt, 


u 


(A.6) 


which, in principle, is much stronger than the one we have in which, we get observability 
estimates on 


z(0,oj; £ 0 ) dIP ( a ;) 


V' 


This difficulty does not arise in the context of exact averaged controllability where, we 

recover the norm of z 0 , and this yields also estimates on z( 0 , uj; z 0 ) for all uj and in particular 
2 

on f n z(0,uj;z o ) dP(uj). 

But the property of null averaged controllability with initial data independent of uj does 
not seem to suffice to derive the same property with initial data that depend on uj. This is 
an interesting issue for further work. 
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A.2 Variational characterization of the controls of minimal norm 

We have shown the existence of the exactly averaged control (resp. null averaged control, 
approximately averaged control), provided that the adjoint system is exactly averaged ob¬ 
servable (resp. null averaged observable, satisfying averaged unique continuation property). 
These results allow concluding whether a system is controllable in an averaged sense. In this 
section, we give variational characterizations of the controls. Such kind of results not only 
derive characterizations of the controls but also serve for computational purposes. 

Theorem A.4 If the system (2.8) is exactly averaged observable, then the exact averaged 
control for the system (2.1) of minimal L 2 ( 0, T; U)-norm is given by 


u(t) = XE(t ) [ B*(u)z(0,u-,zq)¥(u), 
Jn 

where z 0 6 V' minimizes the functional 


J(z 0 ) 


1 

2 


XE{t) B*(u)z(t,uj]z 0 )dF(uj) 


I/O, / z(0,uj;z o )dF(uj] 


vy 


2 

dt — (2/1, zq)v,v> 


(A.7) 


(A.8) 


Remark A.2 The control given by (A.7) is an average of functions of the form B(-)*z(t, •). 
For each sample point u>, B(uj)*z(t,u) can be chosen to be a control. However, generally 
speaking, it does not steer the initial datum yo to the final one y\. 


Proof of Theorem A.f : Define a functional </(■) : V' —>■ M. as follows: 


J(z 0) 


XE(t) / B*(u)z(t,uj]Zo)dF(uj) 


l/o, / 2:(0,o;;2:o)dP(a;; 


vy 


2 

dt — (|/i, Zo)vy 


(A.9) 


It is easy to see that the functional J : V' —» R is continuous and convex. From (2.9), we 
have that J(-) is coercive. Then, we know that </(•) has a unique minimizer z 0 . Let £(•, •) be 
the corresponding solution of the adjoint system. By computing the first variation of J(-), 
it can be seen that 


(yi,zo)vy = [ XE(t)( f B*(uj)z(t,uj;z 0 )du, f B*(uj)z(t,uj; z 0 )d¥(uj)\ dt 
Jo N Jn Jn 1 u 

+ (yo, [ z(0,uj;z o )dF(uj)) ,\/z Q eV'. 

\ Jn / vy 


(A.10) 


From (A. 10), we know that if we choose the control as (A.7), then 


y(T,uj;y 0 )dF(uj) = y x . 
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Now we prove that u(-) given by (A.7) is the control with minimal L 2 (0, T; f/)-norm, 
which drives the mathematical expectation of the solution of the system (2.1) from y 0 to y\. 
Let us choose z 0 = z 0 in (A.10). We have that 


(yi, Zq) vy 



*(0, u; Zo)dF(uj) 




in 


B*(u>)z(t,u>] £ 0 )dP(w) 



(A.ll) 

Let u(-) be another control which steers the mathematical expectation of the solution to the 
system (2.1) from y 0 to y\. Then we obtain that 


(yi,Zo)v,v = f XE(t)(u(t), f B*(u)z(t,u-,Zo)dP(u )^ dt 
Jo J u 

+(yt>, f z(0,w;zo)dP(«))^, Vz 0 e H. 


(A.12) 


From (A.ll) and (A.12), we get that 


XE(t) / 5*M£(0,w;£o)dP(w) 


dt 


u 


XE(t)(u(t), / B*(u)z(0 : uj; z 0 )dP(uj)) dt 


u 


< 


Xe / B*(u)z(0 : uj;zo)dP(uj) 


L 2 (0,T;U) 


\xeu\l 2 (o,t-u)i 


which implies that 


Xe / B*(u)z(0,uj]Zo)dF(uj) 


L 2 (0,T-,U) 


< \Xeu\l 2 {0, T -U)- 


□ 


Remark A.3 The control is the solution of the adjoint system with the final datum which 
is the minimizer of the quadratic, convex and coercive functional J(-) in V'. One can use 
numerical methods such as implementing gradient like iterative algorithms to compute it (see 
[15] for example). However, one will meet the same difficulty as employing this method to 
compute the control for the exact control problems of PDEs. 


Similar to Theorem A.4, we can prove the following result. 

Theorem A.5 If the system (2.8) is null averaged observable, then the null averaged control 
of minimal L 2 (0, T ; U)-norm is given by 


u(t) = XE(t) / B - (w)z(0,u>;2 0 )<fflV), 


(A-13) 


where z 0 e V' minimizes the functional 


j ( z o) = t; [ XE(t) I B*(uj)z(t,uj;zo)dP(uj) dt + (y 0 , f z(0,uj; z 0 )d¥(uj)) 

z Jo Jn u N Jn 1 y y 


(A. 14) 
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Theorem A .6 Suppose that the system (2.8) satisfies the averaged unique continuation 
property. Then for any e > 0, the approximately averaged control is given by 


u(t) = XE(t) / B*(u)z(t,u)dP(u). 


(A.15) 


Here z is the solution to the adjoint system (2.8) corresponding to the datum z o G V' which 
minimizes the functional 


Je(zo) — 


1 


XE(t) / B*(uj)z(t,uj;zo)dF( y uj) dt - (y u z 0 )v,v 

J o U 


Hoi I z(0,u-,z o )dF(uj))^ / + e\z 0 \ V o 


(A.16) 


Borrowing some idea in [54], we can prove the following result stronger than Theorem 
A.6. 

Theorem A.7 Suppose that the system (2.8) satisfies the averaged unique continuation 
property. For any £ > 0 and any finite dimensional space X C V, the solution to (2.1) 
with the control 

u(t) = XE(t) [ B*(u)z(t,uj)dF(u) (A.17) 


satisfies that 


Vi- y(T,u;y 0 )dF(uj) < e, H x yi = ITv / y(T, w; y 0 )d¥(uj). 
Jo V Jo 


(A.18) 


Here U x denotes the orthogonal projection from V to X and z(-) is the solution to the adjoint 
system (2.8) corresponding to the final datum z 0 G V' which minimizes the functional 


= - j XE(t) B*(u)z(t,uj;z 0 )dF(uj) dt - {y u z 0 ) V y 


u 


2 /o, [ z(0,u]Zo)dP(u)\ +e\(I -n* x )z 0 \ V ’- 

Jn 1 vy 


(A.19) 


Proof: Clearly, J E (-) is continuous and convex. We only need to show that it is coercive. 
For this, we prove that 

J £ (z o) —> oo as |zo|v' —*■ oo. (A.20) 

Let {zl)f =1 C V' be a sequence such that \zq\v —> oo as j —* oo. Put z 3 0 = | z 3 0 \ v )z 3 0 for 
j G N. Then 


z o 

1 

1 

\ z o\V' 

2 


2 

C 

U 


+ (yo, / z(0,uj-z 3 o )dF(uj )) +e\(I -Il* x )z 3 0 \ v/ . 

Jn 1 vy 


(A.21) 


If 


lim / XE(t) / B(u)*z(t,u; z 3 0 )dP(u) 
j—>co J 0 Jfl 
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dt > 0, 


u 













then we get from (A.21) that 


J £ {z J 0 ) —* oo as j —> oo. 
Hence, we only need to consider the case that 


lim / Xs(t) / £(w)*2(t,w;3S)<flP(w) 
j-s-ooJo Jq 


dt = 0. 


u 


Since {^olyLi is bounded in V', we can find a subsequence of it, which is still denoted by 
{Aol^Li if there is no confusion, such that z? Q converges to some zo in V' weakly. Thus, we 
have that xe(-)/qB(uj)*z(-, uj; z J 0 )dF(uj) converges to xe(-) J n B(co)*z(-,co] zo)dP(uj) weakly in 
L 2 (0, T; U). Then, we have that 


XE(t) / B(u)*z(t,u-,zo)dF(u) dt < lim / xe{-) / B(uj)*z(t,uj] z 0 )d¥(uj) 


u 


3 ^oo J o 


dt = 0, 


u 


which implies that Xe J n B(lo)* z(t,co] z 0 )dF(co) = 0 in L 2 (0,T;?7). Thus, we know that 
zq = 0. Since X is finite dimensional, we have that 


lim |(J - n* x )zl\ v > = 1. 

J->-oo 


Therefore, 


U4) 


> lim ( - (yi,z 3 0 ) V y' + (y 0 , / z(0,u] z 3 0 )dF(u)) +e) = e, 

j —^oo ' ' J Q ' V,V' / 


lim 

J->oo \Zq\v j^o o 

which implies that J £ (z 3 q) —> oo as \zq\v' ~> oo. By this, we get the coercivity of J e (-). Hence, 
we know that there is a minimizer zq of J £ (-). 

For any 5 > 0 and zo E V', we have that 

0 < ^ [ J £ (z 0 + Szo) - J £ (z 0 )\, 

which implies that 


e\(I -U*x)z 0 \v> < / XE(t)( / B(uj)*z(t,u})z 0 )dF(uj), / B(u)*z(t,uj; z 0 )dF(uj)) dt 


u 


-{vu Zq)v,v> + (yo, / z(0,uj-,z o )dF(uj) 


vy 


Similarly, we have that 


e\(I -n* x )z 0 \ V ' > / XE(t){ / B(u})*z(t,cu;z 0 )dF(uj) : / B(uj)*z(t,uj; z 0 )dF(uj)) dt 


u 


— < 2 / i , zq)v,v + (yo, / z(0,uj-,z o )dF(uj) 
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vy 









Hence, for any zq G V 7 , 


e\(I ~ n^)^o|v' > [ XE(t)( [ B(uj)*z(t,uj;z 0 )dF(uj), [ B(u})*z(t,u; z 0 )dF(u})\ dt 

.In ' ./o ./o ' U 


~(yi,Zo)v,v' +{Vo, / ^(0,a;;^o)^IP(w) 


V,V' 


This, together with (A.2), implies that for any z 0 G V', 

(*o,yi- [ y(T,u;y 0 )) < s\(I^U* x )z 0 \ 

N Jn 1 vy 

Hence, we get (A. 18). 


V'- 
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